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Optical Heterodyne Experiments with 
Enclosed Transmission Paths 


By O. E. DeLANGE and A. F. DIETRICH 
(Manuscript received September 13, 1967) 


Experiments to determine the feasibility of using optical heterodyne 
recewers at the end of long and complicated transmission paths showed 
that nearly ideal performance is possible. The paths were enclosed, thus 
avoiding atmospheric turbulence. We set up two receivers, one coherent 
and the other noncoherent, and compared their signal-to-noise per- 
formances and stability. The degree of beam distortion resulting from 
transmission over a path was determined by comparing the path losses as 
measured by the two receivers. Heterodyne detection was nearly 80 per cent 
efficient for a four-mile path involving 63 reflections and transmission 
through a considerable number of other optical components. Stable hetero- 
dyne operation was observed. 


I. INTRODUCTION 


It is well known that efficient operation of a light heterodyne re- 
ceiver requires a high degree of phase coherence between the received 
signal beam and the local oscillator beam. For example, if the angle 
between the received beam and the local beam is in error by as little 
as one second of are there can be considerable loss of signal. Also it 
is obvious that if the beam were translated by a beam width the de- 
tection efficiency would be reduced to zero. These requirements are so 
critical that many declared successful operation of a heterodyne de- 
tector at the end of a long transmission line to be very doubtful. 

Experiments on the transmission of light beams through the open 
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atmosphere have shown that atmospheric turbulence might render the 
beam unsuitable for heterodyne detection.2, However, the medium 
which we visualize for long-distance optical transmission consists of 
a series of lenses or mirrors for guiding the beam through an enclosure 
which isolates it from its environment. It is quite likely that success- 
ful guided transmission will require some type of beam position 
servo control in some of the guiding elements. If we provide this con- 
trol, the position of the received beam can be stabilized so that an 
efficient heterodyning adjustment, once made, will be maintained. 
Stabilizing the position of the received beam automatically stabilizes 
its angle of arrival. During some previous experiments’? we observed 
that single-mode transmission appeared to be preserved even after 
many reflections. Because of these favorable considerations, we de- 
cided to set up an experimental heterodyne receiver at the end of the 
longest transmission path available and observe its operation. 


Il. THE EXPERIMENTS 


2.1 Transmission Path 


The path consisted of a horizontal aluminum tube six inches in 
diameter and 100 meters long. We put eight concave mirrors, one 
inch in diameter, at one end of this pipe and seven at the opposite 
end to transmit a laser beam for eight round trips through the line. 
The total distance was one mile through enclosed atmosphere. 


2.2 Hxperimental Arrangement 


Fig. 1 shows the experimental arrangement in block form. The laser, 
which operated at 0.634, was in a cavity only 20 cm long. This means 
it could oscillate at only one frequency as long as it was kept tuned 
near the center of the doppler line. Temperature stabilization of the 
cavity was sufficient to assure single-frequency operation without au- 
tomatic frequency control. Power output was approximately one 
milliwatt. 

The laser ouput was amplitude-modulated a few percent by apply- 
ing an 84 KHz sine wave to a modulator which consisted of a suitably- 
oriented rod of KDP followed by an analyzer. The combination of a 
lens and a concave mirror constituted a telescope for adjusting the 
beam diameter and radius of curvature properly for launching it 
into the pipe line. The flat mirror between the lens and the curved 
mirror served only to fold the path back on itself. See Fig. 1. The tilt 
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Fig. 1 — Experimental arrangement, one-mile path. 


of the flat mirror was controlled by voltages applied to two piezoelec- 
tric cylinders. This was part of a servo system used to assure that the 
beam was properly directed to arrive at the center of the mirror at 
the far end of the line. 

After eight round trips through the line the emerging beam, whose 
position was also servo-controlled, was directed along two separate 
paths by a beam splitter. One path included a photomultiplier func- 
tioning as a direct detector for recovering the 84 KHz modulation. 
The second included a telescope for reducing the beam diameter and 
converting it to a nearly plane wave. This wave was combined with 
the beam from the local laser on a second beam splitter. One output 
from this beam splitter was focused onto a photodiode detector. A 70 
MHz intermediate frequency was obtained from this detector as a 
result of the beating of the two signals. The amplified output of the 
70 MHz receiver was applied to a diode detector from which the 84 
KHz modulation was recovered. A 70 MHz frequency discriminator 
furnished an error voltage for automatic frequency control of the 
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local oscillator laser. The amplified error voltage was fed back to a 
piezoelectric mirror mount in such a way as to maintain the frequency 
difference between the two lasers at 70 MHz. 

The beam emerging from the line was centered between four sen- 
sors as shown on Fig. 2..The error voltage from these sensors was 
amplified and applied to piezoelectric cylinders which controlled the 
tilt of the last mirror at the far end of the line and thus kept the 
received beam centered. The remaining mirrors in the line were not 
servo controlled but had manual electrical tilt adjustments. 

The direct detector consisted of an RCA 7326 photo-multiplier 
tube with a circuit in its anode tuned to pass 84 KHz. Because of 
the gain in the multiplier section of this tube it was possible to make 
the signal and shot noise power so much greater than thermal noise 
that the thermal noise could be neglected. The cathode of this tube 
is large enough to make signal output largely independent of beam 
position for moderate amounts of beam translation. The heterodyne 
detector was a Texas Instrument LSX 900 photodiode. A transformer 
tuned to a center frequency or 70 MHz coupled this diode to the 
input of a low-noise preamplifier. The transformer and amplifier had 
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Fig. 2— Servo-controlled beam positioners. 
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a bandwidth of approximately 20 MHz. The active area of the diode 
was only 0.01 inches in diameter which meant that the beam had to 
be very accurately centered on it for maximum response. 


2.3 Determining Heterodyne Efficiency 


The basic plan for the experiment was to measure the loss intro- 
duced by the line as indicated by the heterodyne receiver and to 
compare this with the line loss indicated by the direct detector. If 
the transmission loss indicated by the heterodyne receiver is greater 
than the loss measured by a direct detector the difference can be 
attributed to a decrease of detection efficiency for the heterodyne re- 
ceiver resulting from distortion of the wavefront of the beam by 
transmission through the line. Such distortions should not affect the 
direct detector. 

This can be stated differently by noticing that transmission over 
the long path results in some of the light power being converted into 
higher order modes by any existing line imperfections. That power 
converted to very high order modes cannot propagate in the line and 
is lost. Some of the lower order modes can propagate and will be 
detected by the direct detectors, whereas the heterodyne detector 
can respond to only one mode so that any power not in that mode 
is lost as far as this detector is concerned. 

If transmission of the beam through the required external optical 
components produced wavefront distortion before the beam was 
launched into the line, the distortions introduced by the line might 
be masked. Because of this we need to determine the absolute detec- 
tion efficiency (at least with respect to the beam as it is launched into 
the line). Detection efficiency is defined as the ratio of actual, or meas- 
ured, beat-note current to the current calculated for an ideal detector. 
This efficiency can be measured in several ways. First if the signal 
power, P,, the quantum efficiency, y, and the bandwidth B are known, 
we can calculate the expected carrier-to-noise ratio from the relation- 
ship C/N = 7P;/hfB where h is Planck’s constant and f the light fre- 
quency. By comparing the measured C/N to the calculated value we 
can determine detection efficiency. 

The C/N for the heterodyne receiver was measured by means of 
a detector at the IF amplifier output. With a light signal into the 
receiver the IF gain was adjusted to produce some value of detected 
current to serve as a reference. The signal was then removed and 
the IF gain increased until the detector produced the same amount 
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of current, this time resulting from the noise. This change of gain 
was equal to the C/N. Because the equation applied only to shot 
noise, it was necessary to correct for the presence of thermal noise. 
With the heterodyne gain obtained for this experiment shot noise 
and thermal noise powers were comparable. 

The second method of determining efficiency was to compare the 
amplitude of a low-frequency beat note (produced by beating the 
unmodulated signal with the local oscillator) with the de output of 
the diode detector. It can be shown that, ideally, when the oscillator 
power, P,, is much greater than the light signal power, P,, the ratio 
Ip/Ipg = 2(P;/P.)%. Ip is the peak current at the beat note fre- 
quency. The ratio of measured Ip to calculated I; is the detection 
efficiency. 

The beat note method consistently indicated a higher efficiency 
than the C/N method. Unaccounted-for noise or the cumulative ef- 
fect of errors in determining the light-signal power, quantum effi- 
- ciency, and bandwidth, could account for the pessimistic C/N meas- 
urements. 

We also made some optical tests to determine the heterodyne ef- 
ficiency. Because heterodyne deficiencies result from lack of phase 
coherence between the signal and local oscillator beams, it is evident 
that restricting the area of the superimposed beams to a small region 
near their centers improves coherence. And if this area is made suf- 
ficiently small, perfect coherence must result. Hence, if the combined 
beams are passed through a small iris before being applied to the 
detector, a reference value corresponding to 100 per cent efficiency 
should be obtained. If perfect coherence exists over the total area 
of the beams, increasing the diameter of the iris will not affect the 
efficiency and the manner in which the IF output increases with iris 
area is predictable. If perfect coherence does not exist, increasing 
the iris diameter will reduce the detection efficiency and the IF output 
from the detector will not increase as much as expected. 

For our experiment we used detected direct current as a criterion 
of iris area. Direct current from the detector increases linearly with 
combined local oscillator and signal power. The intermediate fre- 
quency power is proportional to the product of the two light powers 
and should vary as the square of the direct current. Fig. 3 is a plot 
of the results of varying the iris diameter in such an experiment. The 
dashed line shows the ideal variation of IF power with detected direct 
current, that is, the variation which should be expected for perfect 
phase coherence between the two beams. 
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Fig. 3 — Phase coherence tests comparing 2 and 5 dB offsets. 


The solid curve with dots shows actual performance after the 
detector had been adjusted to produce maximum IF power. The per- 
formance is close to ideal. Next, the received beam was tilted from 
the local oscillator beam to reduce IF Power by 5 dB with the iris 
fully open. The curve with circles shows how IF power varied with 
diaphragm opening. The effect of lack of coherence is quite evident 
from the saturation of IF carrier as the iris opening was increased. 
This test was repeated with a 2 dB offset, corresponding to 80 per 
cent efficiency. The curve with squares indicates the results. 


Ill. EXPERIMENTAL RESULTS 


3.1 One Round Trip 

To evaluate the effect of components outside the line, we first 
measured the detection efficiency for only one round trip through the 
line. For this arrangement the signal beam traversed one KDP modu- 
lator, one Glan-Thompson prism, one short focal length lens, two 
spherical mirrors, two flat mirrors, two quarter-wave plates and two 
beam splitters. All of these components were stock items with no 
attempt to obtain exceptional quality. The C/N method indicated 
a detection efficiency of 81 per cent, the beat note method 98 per 
cent; the true value is probably somewhat between. Evidently there 
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was no serious deterioration of performance even though the beam 
returned from the line was noticeably elliptical in contrast with the 
circular beam from the beating oscillator. 


3.2 One-Mile Line-Losses 


By adding 14 more spherical mirrors to the optical path mentioned 
above the path length was extended to eight round trips for a total dis- 
tance of one mile. The efficiency, determined from a number of meas- 
urements of C/N, averaged 81 per cent; the beat note method aver- 
aged 87 per cent. These data indicate that addition of the seven 
round trips to the optical path did not appreciably reduce the detec- 
tion efficiency. This is borne out by the results of another series of 
experiments in which the line losses, as indicated by the heterodyne 
receiver, were found to be the same as the losses measured with the 
direct detector. Table I summarizes the results of this experiment. 

Optical tests of coherence, as we described, were made by passing 
the combined signal and local oscillator beams through an adjustable 
iris. Fig. 4 compares the coherence performance for one and for eight 
round trips through the line. The results agree with those of the other 
tests in indicating no measurable loss of coherence resulting from 
the additional 14 reflections. Within the accuracy of measurement it 
appears that the external optical components produce more beam 
distortion than the one-mile line. 


3.3 Four-Mile Line Configuration 

The round trips through the line were increased by two optical 
circulators to 32 (63 reflections) for a total length of four miles. 
Figure 5 shows the experimental setup. Light from the signal laser 
was first circularly polarized by a quarter-wave plate and then trans- 
mitted through a KDP rod where it was modulated at 84 KHz. One 
component of the modulator output was transmitted through the 
first analyzer prism, the other absorbed. 


TABLE I —One-Mile-Line 


Path Indicated Detection Indicated line Loss (dB) 
(round trips) Efficiency (%) 
Variable Direct 
C/N Beat Tris Heterodyne Det. 
1 81 98 100 — — 
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Fig. 4 — Phase coherence test comparing one and eight round trips. 


The field of the Faraday-rotation isolator was adjusted to produce 
a 45° rotation of the plane of polarization of the amplitude-modu- 
lated output of the prism. A second analyzer prism, at the output of 
the isolator, was set to pass this polarization into the line. A quarter- 
wave plate in the line converted the input to circular polarization. 
After traversing the line for eight round trips the beam was reflected 
directly back on itself and returned to the quarter-wave plate. It 
emerged from the plate linearly polarized in a plane normal to that 
at which it first entered the line. 

Analyzer prism 2 now directed this beam out of one of its side 
ports where a flat mirror reflected the beam directly back on itself 
and thereby launched it into the line for a second time. After 16 
more round trips through the line the beam returned again to the 
quarter-wave plate. Upon emerging from the plate this time it had 
been restored to its original polarization and so was transmitted 
directly through prism 2 into the isolator. The additional 45° rotation 
imparted by the isolator caused its plane of polarization to be 90° 
from that of the original transmission. Upon leaving a side arm of 
_analyzer prism 1, the beam was directed to the heterodyne receiver 
by a flat mirror and two beam splitters. One beam from the first 
beam splitter was applied to a photomultiplier direct detector where 
the 84 KHz modulation was recovered. The other was combined 
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Fig. 5 — Experimental arrangements, four-mile path. 


with the output of the local oscillator and applied to a photodiode 
detector. The 70 MHz output of the photodiode was amplified and 
detected by an electronic detector to recover the 84 KHz modulation. 
The arrangement shown in Fig. 5 made it simple to determine the 
losses introduced by the four-mile line. By placing a flat mirror, 
shown dashed, at the focal point of the lens, light was caused to 
traverse the circulators and other components in the usual manner 
but transmission through the line was eliminated. Thus the line loss 
was the difference between the loss when light passed into the line 
and the loss when the line was blocked by the flat mirror. | 
To determine the effects of local optical components, the absolute 
detection efficiency was determined with only these components in the 
transmission path, that is, with the line blocked by the flat mirror. 
The average of a number of C/N measurements indicated a detection 
efficiency of 86 per cent. The beat method indicated 93 per cent. 


3.4 Four-Mile-Line Losses 


Table II summarizes the results of this experiment. The efficiencies 
and line losses that the table shows are averages obtained from a 
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number of measurements. Fig. 6 shows the variable-iris data. Ad- 
ding the line decreased the detection efficiency from 86 to 78 per cent, 
as indicated by the change of C/N. This is consistent with the line 
loss measurements and with the variable-iris check of coherence. The 
two curves of Fig. 6 show some loss of coherence after transmission 
over the longer path. 

From the results of these experiments one might conclude that the 
four-mile line losses for a heterodyne receiver were about one decibel 
greater than for a direct detector. However, our measurements were 
not accurate enough to make this a reliable figure. It is likely to be 
somewhat pessimistic. In any case, it appears that heterodyne detec- 
tion efficiencies near 80 per cent are obtainable even with light signals 
which have been transmitted through this long and complicated opti- 
cal system. 

We used no really elaborate beam-alignment or beam forming pro- 
cedures for these experiments. We made the signal beam approxi- 
mately plane-parallel by means of a telescope. The local oscillator 
beam remained slightly divergent, just as it emerged from the laser. 
Beam positions and angles were simply adjusted to produce maxi- 
mum IF output from the detector. Although most of the optical com- 
ponents were stock items, the mirrors were rated good to 1/50 wave- 
length by the manufacturers, Laboratory Optical Co. 


IV. STABILITY 


Once it is determined possible to operate a light heterodyne re- 
ceiver at the end of a long transmission line efficiently, the next 
question is whether this efficiency can be maintained for long. Our 
experiments indicate that, with reasonable care, it can. 

There are two types of instability which can cause the system to 
malfunction: frequency variation and beam wander resulting from 


TasBiE IT —Four-Mile-Line 


Indicated Detection 


Path efficiency (%) Indicated line Loss (dB) 
Variable Direct 
C/N Beat Tris Heterodyne Det. 
Local Compo- 
nents Only 86 93 100 — — 
Local Compo- 
nents Plus 


4-Mile Line 78 87 7.5 6.6 
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mirror tilt. Rough temperature compensation in the laser cavities 
was sufficient to keep our transmitting oscillator operating near the 
center of its doppler line. Automatic frequency control kept the aver- 
age frequency of the local oscillator differing from that of the receiver — 
by 70 MHz. The control voltage for this AFC came from a frequency 
discriminator operating at the intermediate frequency. Acoustical 
vibrations acting on the two lasers produced rapid deviations of the 
intermediate frequency obtained from the detector. These deviations, 
with peak values of a few megahertz, were passed without loss by 
the IF amplifier which had a bandwidth of 20 MHz. We have since 
practically eliminated these frequency shifts by using a different type 
mirror mount in the laser cavities. 

Because all of our components were mounted on high temperature 
coefficient materials such as aluminum, brass, and steel there were 
some deviations in the position of the transmitted and received beams 
resulting from tilts produced by temperature variations. Our line, 
which was made up of single mirrors, was very much less stable than 
a practical line, which would have lenses or pairs of mirrors. Also, 
our experimental line was probably subjected to greater temperature 
changes than an underground line would be.* In any case it is evi- 
dent that some form of beam position servo control will be necessary 
in a long transmission line. 


4.1 One-Mile Line 
For this experiment both transmitting and receiving equipment 
were mounted on the same steel plate at one end of the line. The steel 
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Fig. 6— Phase coherence test comparing local and four-mile paths. 


OPTICAL HETERODYNE PERFORMANCE 173 


plate was on a wooden table which stood on a concerete floor. The 
presence of active components, such as the lasers, on this plate re- 
sulted in rather large point-to-point differences in temperature. As a 
result, this plate was the most unstable part of the transmission path 
and, if uncontrolled, the beam received at the far end of the path 
wandered about, sometimes going off the receiving mirror. 

To stabilize the projected beam, four light sensors were installed 
behind the first mirror at the far end of the line. These sensors gave 
error signals which were amplified and applied to piezoelectric cylin- 
ders which in turn controlled the tilt of one of the launching mirrors. 
See Fig. 2. This system improved stability considerably. 

A similar servo system was set up to control the position of the 
received beam to assure proper alignment of signal and local oscillator 
beams. The sensors were at the exit port of the line and the error volt- 
ages controlled the tilt of the last mirror at the far end of the line. 
With such an arrangement it should be possible to maintain the posi- 
tion of the received beam within a few thousandths of an inch—a 
fraction of a second of are in our line. The remaining mirrors in the 
line had electrical tilt controls but these were all adjusted manually. 

We connected a number of monitors to the system to study its 
stability. The outputs of these monitors were continuously recorded. 
Some of the quantities recorded were: (t) Vertical and horizontal 
error voltages from the output beam sensors. These voltages indicated 
the position of the received beam. (i) The 84-KHz signal recovered 
by the heterodyne receiver. This indicated any fading of heterodyne 
signal which might have resulted from such factors as beam wander. 
(wt). The 84-KHz signal recovered from the direct detector. By com- 
paring (i) and (i) we can measure any excess loss resulting from 
decrease of detection efficiency in the heterodyne receiver. (iv) Trans- 
mitted signal power. (v) Local oscillator power. (vt) Automatic fre- 
quency control error voltage. 

From the recordings and from direct observation it was evident 
that the heterodyne receiver was noticeably less stable than the direct 
detector. However, the fluctuations of output of the coherent detector 
were usually less than +0.5 dB and occurred at low rates which could 
easily be removed by an AGC circuit. We should emphasize that the 
heterodyne detector diode was only 0.01 inch in diameter whereas 
the direct detector cathode was over an inch in diameter and could 
accommodate more movement of the received beam. Except for the 
small fluctuations just mentioned, the heterodyne receiver operated 
as well as the direct detector. The link operated all day, and at times 
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for several days, without any readjustment and without serious loss 
of signal. Most times when signal loss did occur it was the result of 
detuning. A better AFC circuit should take care of this. 

The charts in Fig. 7 indicate the stability during a typical six-hour 
period. The important consideration here is the comparison between 
the signals recovered by the direct detector and by the heterodyne 
receiver. Both follow the same general trend and both are completely 
free of any serious fades. 

At point A, near the end of the chart, the servo loop which con- 
trolled the position of the final signal beam was opened. As a result, 
the beam position shifted enough to produce a 1 dB drop in the signal 
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Fig. 7 — Signal received over the one-mile path. 
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from the heterodyne receiver. The output of the direct detector changed 
much less. At point B, the servo control had been restored and the 
signals returned to their former levels. 

That the servo control system can be simple is evident from Fig. 
8 which shows schematically the control circuit for one coordinate 
of one mirror. The components consist of two photoconductor cells, 
two vacuum tubes and two piezoelectric cylinders. The power re- 
quirements are 300 volts at a few milliamperes. With the vacuum 
tube amplifier it would also be necessary to supply heater power. 
When field-effect transistors capable of handling 300 V become avail- 
able it should be possible to substitute them for the vacuum tubes. 


4.2 Four-Mile Line 


Because of the signal reflexing necessary for this experiment it was 
not possible to control the beam position by servo systems. In spite 
of this the configuration was quite stable. Fluctuations of signal level 
from the heterodyne receiver were only slightly greater than for the 
one-mile transmission but were considerably more rapid. These fluc- 
tuations, which had periods of several seconds, appeared to result 
from reflection of energy from the line back into the laser. 

The isolator, supplied by C. B. Rubinstein of Bell Laboratories 
provided 40 dB of isolation. But because of resonances at very high 
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Fig. 8— Automatic beam positioning control equipment. 
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Q the light power could build up in the line, and 40 dB apparently was 
not enough. Tests indicated that these fluctuations resulted from very 
small changes of length of the transmission path. When part of the 
path was through open air where the beam could be affected by air 
currents the fluctuations were considerably greater. 

The charts in Fig. 9 show typical performance of the four-mile path 
for eight hours. The output of the heterodyne receiver decreased 
about 1 dB during this time. At the same time the output of the di- 
rect detector increased almost 1 dB as a result of increased power 
from the transmitting laser. This would indicate that the heterodyne 
receiver faded by almost 2 dB during this interval. If servo control 
of beam position had been possible it probably would have prevented 
this signal loss. 

We do not wish to give the impression that the stability and relia- 
bility of our experimental setup did more than approach that required 
for a practical operating system. However, our observations lead us 
to believe that the required stability can be attained at reasonable 
expense by well known methods. To begin with, a practical line would 
have lenses or pairs of mirrors as beam directors which are much less 
sensitive to mechanical displacement than the single mirrors of our 
line. The practical line would also be underground where it would 
be much less subject to temperature variations.t Having the beam- 
director mounts made of low-temperature-coefficient material would 
help a great deal. Temperature control of the director environment 
could add stability but should not be necessary if beam position is 
controlled by servos. 

Temperature control of the transmitting and local-oscillator lasers, 
along with the best AFC circuit should solve the frequency stability 
problem. Laser tube life would need to be extended. 


Vv. CONCLUSIONS 


We found that reduction in received signal caused by loss of co- 
herence which resulted from transmission through a complicated sys- 
tem (including as many as 70 reflections) was small, comparable 
with the measuring errors. The losses in a four-mile transmission path, 
as measured by a heterodyne receiver, were about 1 dB greater than 
the losses in the same path measured by a direct detector. 

We also found that even without complicated controls the system 
was sufficiently stable to maintain near-optimum operation for fairly 
long periods of time. However, our observations indicate that for long- 
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distance systems which must operate for extended periods with great 
reliability it would be necessary to use more sophisticated control 
circuits. Stability, aside from frequency control, is simply a matter 
of proper beam alignment. This alignment is necessary to obtain 
reliable transmission through the line so that the heterodyne receiver 
requirements are little more severe than the requirements for any 
other type of receiver. 

Although most likely there would be some improvement in efficiency 
through using better optical components, the stock components that 
we used gave quite satisfactory performance. 

It is interesting that even though it is visually evident that the 
mirrors we used had deteriorated and accumulated dust over a period 
of years, the line loss measured now is 0.85 dB per mile in com- 
parison with an average of 0.9 dB per mile measured by a different 
method when the mirrors were new. That is, any change of line loss 
is less than the error of measurement. 
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Measured Beam Deformations in a Guide 


Made of Tubular Gas Lenses 


By P. KAISER 
(Manuscript received September 21, 1967) 


We discuss problems connected with the alignment of gas lenses and 
present experimental results showing the effects of off-axis injection of 
a light beam into a beam waveguide (offset a/r S 0.815, r = radius of 
tubular lens) and of the introduction of various radii of curvatures (R = 
500 m) on the shape and position of the transmitted Gaussian beam. 
The theoretical prediction of reduced beam distortions for nonconfocal 
geometry was confirmed experimentally. We demonstrated low losses 
at small bending radi by using 12 lenses as focusing elements within 
a laser cavity in the bent state. 


I. INTRODUCTION 

Tubular gas lenses have possible use as focusing elements in an 
optical beam waveguide)? Theoretical analysis has shown that 
thermal gas lenses are not ideal; beam distortion from the cumula- 
tive effect of lens aberrations? ® and the beam wandering resulting 
from lens misalignments” * limit the number of lenses that can be 
used in a beam waveguide. Beyond this number, corrective devices 
such as redirectors, refocusers, and mode filters must be used.® *° Re- 
directors sense the position of the beam and introduce deflections which 
tend to realign the beam with the guide axis. Refocusers sense the 
increase of beamwidth and tend to focus the beam back to its ideal 
profile. Mode filters reduce the content of unwanted higher order 
modes. 

Previous experiments indicated difficulties with the accurate align- 
ment of the lenses.1t From this experience we worked out a method 
of alignment by which we achieved satisfactory accuracy and beam 
stability. Subsequently we measured the deformation of the beam 
transmitted off axis and around a bend and compared it with theory. 
We compared beam distortions for nonconfocal Jens spacings with 
those of confocal geometry. 
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Il. EXPERIMENTAL SETUP 


Fig. 1 shows an experimental beam waveguide made of 10 gas 
lenses. Each lens was suspended within a 12-meter-long steel U-chan- 
nel, and kept in place and aligned by two sets of screws as Fig. 2 
shows. Because the lenses were not rigidly connected among them- 
selves except through the U-channel, each could be aligned without 
impairing the alignment of previously-mounted lenses. An air-tight 
interconnection between lenses was accomplished by rubber hoses. The 
gas lenses used are described in Ref. 11. The actual lens elements (6.34 
mm I.D., 15 em long) were mounted within a 5.08 em diameter wave- 
guide tubing (4 meter long), which was connected rigidly to a 5.08 cm 
diameter, 0.5 m long spacer. 

Focal length measurements were performed for a single lens as a 
function of gas flow rate F and temperature rise AT between the in- 
put gas and the heated cylindrical wall of the lens. A collimated laser 
beam was focused through the lens. The beam size as a function of the 
distance d from the geometrical center of the lens was measured (Fig. 
3) and graphically extrapolated back into the lens until it intersected 
the half-power beamwidth of the collimated beam. The relative posi- 
tion of the intersection was a function of F and AT: for low flow rates 
the intersection occurred before, for higher flow rates beyond the 
geometrical center. For constant gas flow rate an increase in tempera- 
ture shifted the intersection toward the beginning of the lens. This 
behavior had been predicted by theory. * The position of the in- 
tersection was obtained with limited accuracy because it was not 
possible to follow the converging beam too far back into the lens. 
Hence, the focal length was determined from the distance d between 
the waist and the geometrical center of the lens. The average width 
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of the collimated beam was w and the width of the waist w,. From 
this the focal length was computed according to 


eerea er 
I= T= way 
Fig. 4 shows the focal length as a function of temperature rise for 
different gas flow rates. For F = 1.5 L/min and AT = 103°C, the prin- 
cipal plane coincided with the geometrical center of the lens. We used 
the resulting focal length of 50 em for the confocal operation of the 
beam waveguide. There was no noticeable difference between the 
focal lengths of the vertical and horizontal plane. 

To guarantee the proper matching, we used two gas lenses as the 
focusing elements in the confocal laser cavity. At first the laser 
cavity was mounted on an optical bench, separate from the trans- 
mission line of 10 lenses. However, the change of the target spot 
position in time, probably caused by motion of the respective sup- 
ports, was not tolerable. We achieved much better stability by mount- 
ing the launching system, together with the transmission line, into 
the U-channel itself. 

In order to guarantee optimum alignment of the laser cavity, we 
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Fig. 3— Beamwidth of collimated and focused light beam as a function of 
distance from a tubular gas lens. 


sent a reference beam through the channel on a predetermined path 
and with its help the optical axes of the two gas lenses of the laser 
were aligned. We placed the glass tube of the laser between the lenses 
and adjusted for minimum scattering of the passing reference beam. 
We placed a plane mirror (mirror 2 in Fig. 1) at the far end from 
the reference laser, so that the energy was coupled back into this 
laser and resonance could be observed by an increase in power output. 
The transmittance ¢ of this mirror was 0.95 percent. After mounting 
a plane parallel mirror at the other end of the laser (mirror 1 with 
t = 0.20 percent), oscillations could be obtained easily and the ref- 
erence beam was shut off. As intended, the aim of the assembled laser 
was close to the predetermined path, and undulations of the beam 
within the cavity were kept to a minimum. The first gas lens of the 
transmission line was closed on one side with a thin glass plate 
mounted at the Brewster angle to block the gas flow. It did not dis- 
place the beam noticeably, a fact which was important for the even- 
tual removal of the plate in order to form a single laser cavity out of 
a larger number of gas lenses. 

The accurate alignment of the lens required the knowledge of the 
beam position inside the lens. For a null indicator we used a probe 
consisting of a differential detecting device with 4 photocells arranged 
as quarter segments of a circle as Ref. 11 describes. In this way, we 
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avoided several difficulties such as calibration inaccuracies, dependence 
upon absolute intensity, and relative low sensitivity for small dis- 
placements of the beam from the probe’s center. We performed the 
alignment as follows. We oriented the cold lens in such a way that 
the probe indicated null readings at each end of the tubular lens. 
Four external indicators marked the horizontal and vertical position 
of the lens within +1 mil. With the probe removed, we determined 
the position of the beam on a target 6 to 12 meters from the lens. 
The target was a micropositioner containing a photodiode and a metal 
plate with a 3-mil hole in front of it. The position of the beam maximum 
intensity could be defined to approximately +8 mils. Then we heated 
the lens to normal temperature. This deflected the beam downward, 
since the optical axis is below the geometrical axis of the tube because 
of the gravity effect. By means of the vertical indicators, we introduced 
a parallel displacement of the lens until the target spot coincided with 
the position before heating. By doing this, we made sure that the 
heating process of the lens did not deflect the beam. Henceforth, its 
trajectory coincided with the optical axis. 

Because of symmetry, no horizontal adjustment of the lens should 


140 


120 





AIR FLOW RATE: 
1.0 L/MIN 








80 


f iN cm 


60 





40 








20 














9) 20 40 60 80 100 120 
AT IN DEGREES C 


Fig. 4— Focal length of the tubular gas lens as function of temperature rise 
between wall and inlet air. 


184 THE BELL SYSTEM TECHNICAL JOURNAL, FEBRUARY 1968 


have been necessary. Even though this was true for most of the lenses, 
a substantial horizontal displacement had to be introduced in some 
of them, apparently to compensate for lens inaccuracies. The vertical 
displacement required for the compensation of the gravity effect varied 
between 13 and 23 mils, with values around 17 occurring most. The 
horizontal displacement necessary, if all at, amounted to a few mils 
and only in one case reached 12 mils. 

One particular difficulty arose in carrying out this procedure. By 
inserting the probe, a back pressure was established and the reduced 
air flow resulted in higher temperatures in the preceding lenses. As a 
consequence, the ensuing variation of the optical axes of these lenses 
deflected the beam from its original target spot. We could make tem- 
perature variations negligible by fast probe insertion. We neglected 
the beam variation caused by reduced air flow and its influence on 
the focal length. Air conditioning machinery vibrated the building 
and caused random displacements and amplitude fluctuations in the 
laser beam. This was particularly detrimental for the centering pro- 
cedure and prevented our using a more sensitive range in the probe’s 
differential amplifier. The fluctuations disappeared after the motors 
were stopped for a short period. Taking the various factors into account, 
we estimated the alignment accuracy of the lenses’ optical center to 
be within +3 mils, which is 4 of the theoretical 1/e beam halfwidth 
w(= 17.75 mils) at each lens. 


II. OFF-AXIS INJECTION 


We injected the light beam off axis by passing it through a tilted 
plane parallel glass plate 83.5 mils thick. Figs. 5 through 8 show the 
power profiles of the beam transmitted through the guide on axis and 
with different horizontal and vertical offsets. We took about 25 points 
for each power profile and drew a smooth curve through these points 
as Fig. 5 shows. The target was 6.25 meters from the center of the 
last lens. For a beam injected on axis, the half-power beamwidth 
changed from day to day. It varied between 0.13 and 0.16 inch verti- 
cally, and between 0.15 and 0.17 inch horizontally. The center of the 
beam on the target also moved slowly a few tens of mils as a function 
of time. We found no satisfactory explanation for this behavior. 

We normalized both the vertical and horizontal beamwidths for 
off-axis injection with the horizontal beamwidth of the on-axis beam. 
We obtained this reference width at the beginning of each series of 
measurements. We normalized the distance between the maximum 
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Fig. 6 — Vertical power profile for horizontally-displaced beam. 
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intensities for off-axis and on-axis beams with the radius of the tubu- 
lar lens and set the maximum intensities to unity. As could be ex- 
pected from theoretical considerations, the deformations in the verti- 
cal and horizontal plane were independent of each other to the first 
order. However, we noticed a stronger deformation of the vertical 
beam shape in case of larger horizontal displacements, as compared 
with the horizontal beam shape for large vertical displacements (see 
Figs. 6 and 8). The most important result shown in Figs. 5 and 7 is 
the increase in beamwidth for growing offsets and eventual beam 
deformations for large offsets. Fig. 9 shows the increase in horizontal 
and vertical beamwidths relative to their respective values at zero 
offset. Points obtained for different experimental runs and, in general, 
resulting in different beamwidths, fitted well into the graph, once they 
were related to their respective horizontal zero offset widths. The 
uncertainty of the result owing to limited accuracy of the measure- 
ments and ambiguity of interpolation becomes larger for increasing 
beamwidths and is indicated in the graphs. The growth appears to be 
quadratic to the first order. Whereas this dependence 1s approximately 
equal for left, right, and high injected beams, the growth is larger for 


TUBULAR GAS LENSES 187 


the beam injected low. We found no satisfactory reason for this dis- 
crepancy. 

The good symmetry in the position and the power profile of hori- 
zontally deflected beams appear to be the result of good alignment 
and horizontal symmetry of the lenses. The asymmetric deforma- 
tions of the lenses owing to gravity apparently result in asymmetry 
of the position and profile for vertically-displaced beams. 

For confocal geometry the normalized width w, of a light beam 
which undulated with amplitude a about the axis of a sequence of n 
lenses is given by?® 


2 
w,=1+3 in(4) 
w 


where § is a measure of the focal length distortion. At distance w from 
the center of the lens the relative increase of the focal length is 8. D. 
Marcuse’s analysis of the tubular gas lens‘ yields a value $ of approxi- 
mately 0.012 for the geometry, gas flow rate, and temperature rise 
being considered. The theoretical prediction for the beam size increase 
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Fig. 8 — Horizontal power profile for vertically-displaced beam. 
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according to above equation corresponds closely to the experimental 
results. (See Fig. 9.) Ten per cent beam size increase occurs after 10 
lenses with a beam displacement a/w of only 0.7. 


IV. BENT BEAM WAVEGUIDE 


Because the lenses were mounted within a U channel it was easy 
to put bends of varying radii into the vertical plane. We lifted the 
channel to calculated heights for this purpose. The curvature started 
at the first lens of the optical guide (third lens in Fig. 1). The laser 
cavity with its two gas lenses remained straight. The U-channel was 
sufficiently elastic to permit following the desired curvature. The light 
beam entered the bend on the optical axis of the first guide lens. Fig- 
ures 10 and 11 show the beam profile transmitted through the line for 
five different radii between 2500 and 500 meters. The half-power 
widths are normalized with 0.425 cm. No significant change of the 
Gaussian beam shape occurred down to a radius of 834 meters. A small 
hump appeared on one side of the beam at the 625 meter radius, but 
we measured no noticeable power loss. At 500 meters, however, the 
beam showed characteristic diffraction rings and a power loss of ap- 
proximately 1 dB. 
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Fig. 9 — Relative increase of beam width as a function of horizontal and verti- 
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Fig. 10 — Vertical power profile as a function of radius of curvature. 


According to the theory of a circular bend,” a beam injected off 
axis and with a particular tilt does not undulate, whereas on-axis in- 
jection causes undulations. We found that a small offset (without tilt) 
of approximately 0.35 mm for a 500-meter radius resulted in optimum 
transmission (power loss approximately 0.2 dB). For a parallel dis- 
placement without tilt a theoretical offset of 0.6 mm results in un- 
dulations around the ideal trajectory with minimum amplitude. As a 
consequence of this, the distortions are also minimized. Just as for 
vertical displacements, the horizontal shape did not change signifi- 
cantly for vertical bends. 

In order to confirm the low-loss guidance properties of the 10-lens 
bent beam waveguide, we transformed the 2-lens laser cavity and the 
10 guide lenses into a 12-lens laser, whereby the discharge tube of the 
2-lens laser remained the only gain medium (see Fig. 1). Lasing was 
obtained at R = 834 meters with power output equal to that of the 
laser with only 2 lenses. The alignment procedure was as follows. 
With the 2-lens laser oscillating, we mounted a plane parallel mirror 
behind the 12th lens (mirror 3 in Fig. 1) and adjusted this mirror for 
resonance. Then we removed mirror 2 as well as the Brewster angle 
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Fig. 11 — Horizontal power profile as a function of radius of curvature. 


window between the laser and the guide and connected the two adjacent 
lenses by a hose. The 12-lens laser started lasing immediately. We re- 
placed mirror 3 by mirror 2 in order to obtain identical operating con- 
ditions. Subsequently, the 12-lens laser was deformed into a circle with 
a radius of 625 meters (the first two lenses remained straight). Even 
though we did not obtain stable lasing, the proximity to oscillations 
existed, since small changes of air flow were enough to produce light 
flashes. 

By changing the alignment of the mirrors, we created various 
paths through the 12 lenses of the laser cavity. This affected the field 
configuration of the fundamental mode. We observed strongly deformed 
beam shapes whose single mode configuration could be confirmed by 
placing a circular aperture into the path of the beam. We recorded no 
change of these shapes when the aperture was closed until laser action 
vanished. 


V. BEAM WAVEGUIDE WITH NONCONFOCAL SPACING 


Theoretical work by E. A. J. Marcatili?® indicated that distortions 
of the Gaussian beam passing through a sequence of lenses are reduced 


TUBULAR GAS LENSES 191 


if a nonconfocal geometry is chosen. To check this prediction experi- 
mentally, we aligned a beam waveguide consisting of 10 gas lenses with 
nonconfocal spacing (a distance-to-focal-length ratio (d/f) of 1.82). At 
the same time, we could make several improvements in the alignment 
procedure as a consequence of the experience gathered with the beam 
guide aligned before. We strongly reduced the influence of building vi- 
brations by working farther away from the maintenance machinery at 
the far end of the waveguide corridor. The ensuing reduced fluctuations 
of the light beam allowed a more accurate coaxial alignment of the 
tubular lens with the light beam. Also, we eliminated the negative in- 
fluence of the back pressure due to probe insertion on the alignment by 
leaving a small air gap between the lens (which was to be centered 
around the beam) and the preceding lens. This permitted the gas to 
escape. After the centering procedure we connected the two lenses air 
tight and then compensated for the gravity effect. Because it was a 
different set of gas lenses than in the experiment before, the vertical 
displacement to compensate for the gravity effect varied between 17.5 
and 20.5 mils, but reached extremes of 15 and 26 mils. 

We left the laser cavity untouched, and adjusted the temperatures 
of the guide lenses to yield focal lengths of 60 and 50 em, correspond- 
ing of d/f = 1.67 and 2.0, respectively (Figs. 12 and 13). This did not 
deflect the on-axis beam in the horizontal direction. However, small 
deviations in the vertical position are apparently the result of the tem- 
perature dependence of the gravity effect. The horizontal and vertical 
beam widths of the on-axis beam for confocal geometry were virtually 
identical with those we measured previously. For nonconfocal opera- 
tion, the horizontal beamwidth of the on-axis beam remained approxi- 
mately constant, whereas the vertical width became larger for noncon- 
focal geometry. We recorded good symmetry of the target beam posi- 
tion with respect to symmetric off-axis injection of the beam in all 
three cases. The normalized offsets are the same as in Figs. 5 through 
8. Again we measured the power profile 6.25 meters behind the last 
lens of the guide. However, we chose arbitrary units at this time 
whereby the power into the guide was held approximately constant. 
We could not confirm the asymmetry observed in Fig. 7 for small 
vertical offsets. 

If you compare the power profiles of light beams passing through a 
sequence of lenses with different distance-to-focal-length ratios, you 
notice a marked increase in beam distortions as you approach confocal 
geometry. This is true for beams undulating horizontally as well as 
vertically. This behavior was predicted by theory?® and confirms that 
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in beam waveguides with nonconfocal spacing there is less mode con- 
version than in those with confocal spacing. 


VI. CONCLUSION 


We achieved accurate alignment of a sequence of gas lenses. We 
found that a beam undulating in an optical waveguide increases its 
beamwidth as a function of the initial offset, but maintains its Gaus- 
sian profile up to a certain maximum displacement. With 10 lenses 
we recorded maximum offsets of approximately 1.0 to 1.5 (normalized 
with the theoretical 1/e beamwidth w = 0.45 mm at the positions of 
the lenses). Larger offsets resulted in a change of the Gaussian beam 
shape and the appearance of diffraction rings. Observed symmetry in 
the position of the maximum intensity and the power profile for left 
and right horizontal displacements indicate symmetry of the lenses 
with respect to the vertical plane. Unsymmetrical vertical deformations 
are apparently the consequence of gravity aberrations. 
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Fig. 12— Horizontal power profile for. horizontally-displaced beam (arbitrary 
units, offset normalized as in Fig. 5). (a) focal length: 60 cm (d/f = 1.67). (b) 
focal length: 50 em (confocal). 
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Fig. 18 — Vertical power profile for vertically-displaced beam (arbitrary units, 
offset normalized as in Fig. 7) (a) focal length: 60 em (d/f = 1.67). (b) focal 
length: 50 cm (confocal). 


Thus, averaged aberration properties of gas lenses can be ascertained 
from offset injection of a light beam into a small sequence of lenses. 

We demonstrated effective guidance of a light beam around a circular 
bend. The minimum radius below which changes of the beam shape oc- 
curred was approximately 625 meters for 10 lenses. We confirmed the 
low-loss guidance properties by using 12 lenses as focusing elements 
within a laser cavity in the bent state. 

We experimentally confirmed the theoretical prediction of reduced 
beam distortions for nonconfocal geometry. Thus, nonconfocal spacing 
of the lenses is more advantageous in optical transmission line design 
than confocal spacing. 
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Silicon Schottky Barrier Diode with 
Near-Ideal LV Characteristics 


By M. P. LEPSELTER and S. M. SZE 
(Manuscript received October 19, 1967) 


Metal-semiconductor diodes with near-ideal forward and reverse I-V 
characteristics have been fabricated using PtSi contacts and diffused 
guard rings. Typically, for a device with an area of 2.6 X 10°° cm made 
on an n-type (111) oriented, 0.85 ohm-cm silicon epitaxial substrate, the 
forward current follows the expression I, = I, exp (qV/nkT) over eight 
orders of magnitude in current with I, = 10°" A and n = 1.02. The 
reverse breakdown is sharp and occurs at the theoretical breakdown voltage 
of p'n silicon junctions of the same n-type doping. The premature break- 
down observed in nearly all previous Schottky barrier diodes has been shown 
to be caused by electrode sharp-edge effects. Besides giving sharp breakdown 
voltage, the guard ring also eliminates anomalously high leakage currents, 
yet stall retains the fast recovery time characteristic common to other Schottky 
barriers. Typically, the recovery time measured at 10 ma ts less than 
0.1 ns, the resolution of the measurement. 


I. INTRODUCTION 


Metal-semiconductor (or Schottky) barriers have been studied ex- 
tensively in the past few years. The main concerns are the barrier height 
and the current transport in the metal-semiconductor system; the 
former has been reviewed recently by Mead," and the latter discussed 
by Crowell and Sze.? Most emphasis has been placed on the forward 
current-voltage characteristic, which relates intimately to the elec- 
tronic applications of Schottky diodes such as varistors for logic gates 
and microwave downconverters. The reverse current-voltage charac- 
teristic of a conventional planar Schottky diode,? usually has had 
greater leakage current and lower breakdown voltage than a diffused 
p-n junction. Because of this “soft” reverse characteristic, Schottky 
diodes have not been considered for power application, as well as 
IMPATT oscillators. 
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It should be pointed out, however, that if the reverse characteristic 
can be improved substantially, a Schottky diode with its inherent 
majority transport property, can be used as a high-speed switch in 
which there is virtually no minority carrier storage, and as a high- 
power rectifier or high-frequency oscillator where the region of great- 
est heat dissipation is located right at the metal-semiconductor inter- 
face and therefore the heat can be more readily conducted away. 

The soft breakdown of conventional Schottky diode is not caused 
by the electron tunneling effect but results mainly because of the 
“edge effect”? shown in Fig. 1(a), where high-field concentration gives 
rise to excess leakage current and low breakdown voltage. The tunnel- 
ing effect is ruled out because of the fact that even at very large electric 
fields (~5 xX 10° V/cm), there is negligible contribution of the 
tunneling current component to the total conduction curent,? owing 
to the relatively large electron effective mass in silicon such that the 
tunneling probability of electrons from the metal to the semiconduc- 
tor conduction band is very small. 

This paper presents one method to eliminate this edge effect: the 
diffused guard-ring method. This method has been used on planar 
p-n junction devices‘ to eliminate the junction curvature effect. It is 
shown in this paper that the guard ring improves both the forward and 
reverse charactcristics of the Schottky diode. It is also demonstrated 
for the first time that the measured breakdown voltage of a Schottky 
junction is equal to the theoretical value of a one-sided abrupt p-n 
junction with the same background doping concentration. In addi- 
tion, the surface field effect associated with the metal-inulator- semi- 
conductor structure on the junction breakdown is studied. 


II. EXPERIMENTAL PROCEDURE AND RESULTS 


2.1 Device Fabrication 

In order to study the effect of a junction guard ring on the I-V 
characteristics of a metal-semiconductor barrier, three kinds of struc- 
tures were fabricated on a single wafer as shown in Fig. 1. Fig. 1(a) 
is a planar PtSi-Si Schottky diode alone; Fig. 1(b) is the diffused p-n 
guard ring alone; and Fig. 1(c) contains both a PtSi-Si Schottky 
barrier and a p-n junction which serves as its guard ring. 

Most of the wafers used in the study were n-type silicon materials, 
with 2 to 8 wm thick epitaxial layers, (111) oriented, 0.3 to 1.2 ohm-cm. 
The basic steps of fabricating the Schottky diode with a guard ring 
were: the wafers were first cleaned and degreased; silicon dioxide 
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Fig. 1— (a) Planar PtSi Schottky diode. (b) Diffused p*n junction ring. (c) 
PtSi Schottky barrier with a p*n guard ring as in b. 


layers of about 5000 A were then grown thermally. Diffusion rings 
were then cut in the oxide by the standard photoresist technique. 
After the junctions were boron-diffused, another SiO. layer was 
grown on the surface, and photoresist was used to define the center 
hole for platinum deposition. 

In order to produce a metal-semiconductor interface comparable 
in uniformity to a diffused p-n junction it is required to start with 
a “clean” silicon surface. Even the most efficient chemical methods 
leave several atomic layers of inorgonic films. In this study, the silicon 
surfaces were back-sputtered with argon ions (at a pressure of 20 
microns) using the setup shown in Fig. 2. Cathode 1 was excited with 
an RF oscillator, producing a plasma over the sample, and removing 
any films previously formed. Then a DC voltage was applied to cathode 
2, and a 500 A platinum film deposited. Heating the sample in the 
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Fig. 2— Platinum sputtering system. 


vacuum chamber to 600°C caused the platinum film to react with the 
exposed silicon surface. For a clean surface, the reaction proceeded 
to the most stable phase of the system Pt-Si, forming the compound 
PtSi. This reaction, verified by X-ray analysis, is strongly influenced 
by the actual cleanliness of the interface. 

The platinum layer deposited on the SiO., which did not react with 
the exposed silicon, was then removed by aqua regia. The final overlay 
metal contact was made by evaporation of aluminum and another 
photoresist to isolate each device. The back contact was formed by 
evaporation of gold and alloyed at 400°C. The three different structure 
shown in Fig. 1 were made at the same time on the same wafer by 
using multiple patterns on the photoresist masks. 


2.2 Forward I-V Characteristics 


Figure 3 shows the measured forward J—V characteristics for the 
three structures as shown in Fig. 1. The epitaxial thickness is 2 ym, 
and the background doping is 2.2 x 10°° cm™*. The diffused junction 
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depth is about 1 wm. The guard ring has an area of 2.9 X 107° 
and the Schottky junction of 2.5 X 10°° cm’. One notes that the 
current of the Schottky diode with a guard ring (from now on, we 
call it diode c referring to Fig. 1c) follows the expression J, = J, exp 
(gV/nkT) over eight orders of magnitude with J, & 10° amp and 
n = 1.02. From the saturation current and the area of the diode (2.5 X 
10° cm’) a barrier height of 0.85 eV for the PtSi-Si barrier is obtained, 
in agreement with the previously reported results.° 

For the planar Schottky diode without a guard ring (diode a) it 
is apparent that the J-V characteristic is much inferior and gives 
about two to four orders of magnitude excess current at lower bias. 
The forward J—V characteristic of the p-n junction guard ring is shown 
in Fig. 1(c). It is clear that the guard ring helps to eliminate the edge 
leakage current which has normally existed in a planar Schottky diode. 
By the use of lower resistivity and a thinner epitaxial layer, it is ex- 
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Fig. 3— Forward current-voltage characteristics of devices shown in Fig. 1. 
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pected that one can obtain an exponential current range (with con- 
stant n-value) even larger than the present eight orders of magnitude. 


2.3 Reverse I-V Characteristics 


The Schottky diode with a guard ring also has a superior reverse 
I-V characteristics. For diode a the breakdown voltage is only 5 V 
(at 1 ma). The breakdown voltage for diode b is 27 V. About the 
same breakdown voltage is obtained for diode c which is shown in 
Fig. 4 for two different areas of the Schottky diodes. For small area, 
the space-charge generation and recombination current of the p-n 
junction guard ring dominates the reverse leakage current. For larger 
area, however, the reverse current approaches the ideal Schottky 
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Fig. 4— Reverse current-voltage characteristics of Schottky diodes (type ¢ in 
Fig. 1). For the small diode, the area of the Schottky diode is 2.5 X 10° em? and 
and that of the guard ring is 29 x 10° cm*. For the larger one, the areas are 
2 xX 10° cm? and 1.5 x 10 cm’, respectively. 
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barrier current due to image-force lowering (the dotted line is calcu- 
lated using an image-force dielectric constant of 12).5 The breakdown 
voltages in diodes b and c are in good agreement with the theoretical 
calculation shown in Fig. 5. This figure shows the breakdown voltages 
for Si p*n junctions as a function of the background doping with the 
junction depth as a parameter.® 

For planar diffusion through a photoresist mask, the impurities 
will diffuse both downward into the bulk semiconductor and sideways 
under the oxide layer (see insert of Fig. 5). Thus a junction curva- 
ture is formed near the corner of the diffusion mask. For a one-sided 
abrupt junction, this curvature causes a reduction of the breakdown 
voltage. In the present case with Ng = 2.2 x 10%® cm‘, the break- 
down voltage would be 36 V for a plane junction, and about 25 V for 
a junction with 1 »m curvature,® as shown by the circles in Fig. 5. 
Therefore the breakdown voltages for diodes b and c agree with theo- 
retical expectation for a guard ring with 1 »m junction depth. 

The above result confirms that a Schottky diode with a junction 
guard ring can have a breakdown voltage many times larger than 
one without a guard ring. In this example, breakdown occurred at 
the guard ring; but this result does not imply that the breakdown 
voltage of a Schottky diode is always limited by the curvature effect 
of the junction guard ring. 
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_Fig. 5— Theoretical results for silicon p*n junction with junction depth r; 
(insert) as the parameter. (See Ref. 6.) 
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By proper control of the diffusion process, one can form a linearly- 
graded junction with Vz which is virtually independent of the junction 
curvature,* or one can form a composite junction (see insert of Fig. 6) ; 
that is, the space charge terminates in a graded region on one side of 
the junction and in a uniformly doped region on the other. The com- 
posite junction has a breakdown voltage always larger than an abrupt 
junction with the same junction curvature (7r;) and the same back- 
ground doping (Nz). Some theoretical results of the composite junc- 
tion are shown in Fig. 6. 

In order to make a guard ring which has higher breakdown voltage 
than that of the Schottky diode, the following experiment was per- 
formed. An epitaxial silicon wafer with doping about 2 X 10'° em™ 
and a thickness of 8 wm was used to form the three structures as shown 
in Fig. 1. The diffused guard ring has a doping profile which can be 
approximated by a composite junction with a junction depth of 3.5 um 
and an impurity gradient of about 2.5 X 10°’ em™*. The theoretical 
breakdown voltage for the guard ring is about 50 V and the breakdown 
voltage for the Schottky barrier is about 42 + 2 V (shown in I'ig. 6). 
Thus the Schottky barrier diode is expected to break down first. Figure 7 
shows the I-V characteristics of the Schottky diode with guard ring 
(diode c) which has a forward characteristic similar to that shown 
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Fig. 6 — Theoretical results for silicon linearly-graded composite junction (in- 
sert) with junction depth as the parameter. 
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_Fig. 7— Oscilloscope current-voltage display of the characteristics of a Schottky 
diode with the ideal breakdown voltage (Nv = 2.0 x 10% cm™%). 


in Fig. 3, and a reverse breakdown at about 41 V in good agreement 
with the expected result. 

Figure 8 shows the forward and reverse J—V curves for the three 
structures. For diode a the forward and reverse characteristics are poor. 
For the p-n junction (guard ring alone, see Fig. 8b) the forward knee 
occurred at about 0.5 V and the reverse breakdown at 49 V. For diode c, 
the Schottky diode breaks down first at about 48 V. Because of the 
small junction area there is a large space-charge resistance, R,, , under 
avalanche conditions. The value of R,, is given by’ R,, & [W’/A (2eV.)] 
where W is the width of depletion layer at breakdown, A the junction 
area, e the permittivity, and V,z the limiting drift velocity. In the 
present case with W & 3 wm, A = 2.5 X 107° cm’, and V, & 10’ cm 
per second, the value of R,, is 1.5 KQ, which is in good agreement 
with the slope shown in Fig. 7(c). (The junction temperature effect 
and the series resistance effect are calculated to be small compared 
with the above space-charge effect.) This further confirms that we 
have uniform breakdown of the Schottky diode. When the voltage 
reaches about 50 V, the p-n junction guard ring begins to break down 
similar to case b (Fig. 8b). The space-charge resistance in this region 
is considerably lower because of the larger area of the guard ring. 


2.4 Surface Field Effect Study 

Notice in Fig. 1(c) that in addition to diffused p* guard ring, there 
is a metal-oxide-semiconductor structure surrounding the device. It 
is known that when a reverse bias is applied to the metal electrode, 
the surface space-charge region under the oxide is in nonequilibrium 
condition.® 

It has been shown that a surface field has a profound effect on the 
breakdown of planar p-n junctions.® In order to study this effect on 
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Fig. 8—Oscilloscope display of the current-voltage characteristics for the 
structures shown in Fig. 1. 


the breakdown voltage of Schottky barrier diodes, a separate metal 
electrode on the oxide (that is, the gate electrode) is fabricated as 
shown in Fig. 9. The fabrication procedures were the same as de- 
scribed previously. The samples used were (111) oriented, 1.0 ohm-cm, 
n-type silicon wafers. When a negative gate bias is applied, the sur- 
face field tends to smooth out the field concontration near the junction 
edge (see Fig. 10). Thus the radius of curvature, 7;, is effectively in- 





ae 9—Schottky barrier diode with separate electrodes on PtSi and on the 
oxide. 
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Fig. 10 — Electric field distribution for negative gate voltage. 


creased. This in turn increases the breakdown voltage (the relation- 
ship between Vz and r; has been shown in Fig. 5). 

Figure 11 shows the field profile near the junction as a positive gate 
bias is applied, where the radius of curvature, 1; , is effectively reduced 
resulting in a decreased breakdown voltage. Figure 12 shows the meas- 
ured reverse J—V characteristics as a function of the gate voltage for 
a Schottky diode with a guard ring of 0.4 um junctions depth. As ex- 
pected, the gate voltage does have a profound effect on the junction 
breakdown voltage. Figure 13 shows the measured breakdown voltage 
versus gate voltage for four different junction depths. Notice that the 
breakdown voltages all approach the theoretical value (100 V) as 
—Vg increases. Also at zero bias, the breakdown voltage decreases as 
the junction depth 7; decreases. 

Similar effects are observed on a Schottky diode that has no p-n 
junction guard ring but has a second metal oxide semiconductor over- 
lay as shown on the Fig. 14 insert, where ZrO2 is formed near the 
periphery of the metal, and separate voltages are applied to the 
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Fig. 11 — Electric field distribution for positive gate voltage. 
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Fig. 12— Reverse current-voltage characteristics as a function of the gate 
voltage on a Schottky diode with a p-n junction guard ring. 
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Fig. 13— Measured breakdown voltage (at 1 ma) versus gate voltage for 


Schottky diode with p-n junction guard ring. The junction depths are 0.4, 1.2, 3, 
and 5 wm. 
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Fig. 14— Measured breakdown voltage versus gate voltage for Schottky diode 
with a metal oxide semiconductor overlay near the periphery (insert). 


diode and the metal oxide semiconductor guard ring electrodes. As 
the gate voltage increases in the negative direction, the edge field is 
- gradually reduced, resulting in an increased breakdown voltage (here 
defined as voltage drawing 1 ma current). When Vg increases in the 
positive direction, however, the enhanced edge field causes drastic 
lowering of the junction breakdown voltage as Fig. 14 shows. 


III. SUMMARY 

It has been demonstrated that the breakdown voltage of a Schottky 
barrier diode can be made to approach the theoretical value of one- 
sided abrupt p-n junction with the same background doping. The much 
lower breakdown voltages usually obtained result from “edge effect” 
which can be modified or eliminated by proper use of either a dif- 
fused guard ring or a surrounding electrode. It has also been shown 
that the surface field can have profound influence on the guard ring 
breakdown characteristics. Since any planar diffused guard ring has 
a finite junction curvature (r;), to eliminate the junction curvature 
effect one has to control the diffusion depth and diffusion profile 
such that the breakdown voltage of the guard ring is larger than that 
of the Schottky junction. 


IV. ACKNOWLEDGMENT 


The authors wish to thank D. Kahng for his stimulating discus- 
sions. They also wish to thank R. W. MacDonald for his assistance 


208 THE BELL SYSTEM TECHNICAL JOURNAL, FEBRUARY 1968 


in processing the samples, Mrs. M. H. Read for her X-ray analysis 
of the PtSi films, and A. Goetzberger and R. M. Ryder for helpful 
suggestions. 


REFERENCES 


1 


. Mead, C. A., Metal-Semiconductor Surface Barriers, Solid State Electron, 


9, 1966, pp. 1023-1034. 


2. Crowell, C. R. and Sze, S. M., Current Transport in Metal-Semiconductor 
3. 
4 


Barriers, Solid State Electron, 9, 1966, pp. 1035-1048. ; 
Kahng, D. and Lepselter, M. P., Planar Epitaxial Silicon Schottky Barrier 
Diodes, B.S.T.J., 44, 1965, pp. 1525-1528. 


. Goetzberger, A., MacDonald, B., Haitz, R. H., and Scarlett, R. H., Ava- 


lanche Effects in Silicon pn Junctions II. Structurally Perfect Junctions, 
J. Appl. Phys., 34, 1963, pp. 1591-1600. 


. Sze, S. M., Crowell, C. R., and Kahng, D., Photoelectric Determination of 


the Image Force Dielectric Constant for Hot Electrons in Schottky Bar- 
riers, J. Appl. Phys., 35, 1964, pp. 2534-2536. 

Sze, 8S. M. and Gibbons, G., Effect of Junction Curvature on Breakdown 
Voltage in Semiconductors, Solid State Electron, 9, 1966, pp. 831-845. 


. Sze, S. M. and Shockley, W., Unit-Cube Expression for Space-Charge Resis- 


tance, B.S.T.J., 46, 1967, pp. 837-842. 


. Grove, A. S and Fitzgerald, D J., Surface Effects on pn Junctions: Char- 


acteristics of Surface Space-Charge Region Under Non-Equilibrium Con- 
ditions, Solid State Electron, 9, 1966, pp. 783-806. 


. Grove, A. §., Leistiko, O., and Hooper, W. W., Effects of Surface Fields on 


the Breakdown Voltages of Planar Silicon pn Junctions, IEEE Trans. 
Electron Devices, ED-18, 1967, pp. 157-162. 


Design of a Simulator for Investigating 
Organic Synchronization Systems 


By R. H. BOSWORTH, F. W. KAMMERER, D. E. ROWLINSON 
and J. V. SCATTAGLIA 


(Manuscript received October 24, 1967) 


This paper describes circuits and devices designed for a real time sim- 
ulator that was used to study a network synchronization scheme. The 
simulator represents a system of four geographically-separated stations 
with compound interconnections. The system clock rate is a function of 
the states of each local clock and the interconnecting paths. Component parts 
for the simulation system were built to specifications exceeding require- 
ments of an actual system so that the operating states could be precisely 
determined. The simulator performed reliably during the series of expert- 
ments which provided the data described in a companion paper.” 


I. INTRODUCTION 


The problem of maintaining synchronism in a pulse code modula- 
tion network has been of interest for some time. 

One technique for synchronizing the geographically separated 
clocks of a network is known as organic synchronization. In this 
scheme, the frequency of each clock is controlled to maintain a phase 
lock between its signal and the average phase of signals arriving from 
other stations. 

The next article in this issue describes a design of such controls 
and also describes simulations that demonstrate the practicality of 
the technique. This paper describes the analog simulator used in 
that work. ’ 

The organic synchronizing method compensates the slow phase 
drift caused by discrepancies in the clock frequencies and the phase 
changes brought about by changes of transmission delays. The fast 
phase jitter introduced by changes in the PCM word patterns is not 
the concern of this study.? Therefore, in the simulator, it is convenient 
to use the cycles of a constant amplitude sinusoid to represent the 
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PCM frames. The sinusoids are generated by independent oscillators, 
one situated at each station, and the transmission paths between them 
are simulated by delay lines. 

The particular network simulated is shown in Fig. 1. It comprises 
four stations interconnected with 12 links. Each link has a large fixed 
delay representing a transmission path typical of transcontinental 
distances. A small variable delay is added for simulation of slow 
delay changes. 


II. A STATION DIAGRAM 


The control functions for a single station are illustrated in Fig. 2. 
This station is designated station 7. The other stations are labeled j, 
k, and Ll. In Fig. 2, only the link from station j is shown in detail, but 
similar equipment is provided in all links arriving at each of the four 
stations. The function of the system blocks will be described while 
tracing a signal through the network. 

The sinusoidal output of oscillator j7 comes to station 7 through a 
cascade of fixed delay elements selected to represent a particular path 
length. The last delay element in the path is servodriven and is con- 
tinuously variable under control of an external signal to simulate 
transmission delay variation. 
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Fig. 1— A network for simulating organic synchronization of four geographi- 
cally-separated stations. 
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Fig. 2— Control functions for one of four stations. The illustration shows 
station 1 and the simulated transmission links connecting with station j. 
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At station 7 the signal is fed to a second continuously variable delay, 
the compensating buffer store, that is servodriven to maintain the 
required phase lock. The drive for this delay is derived from the prin- 
cipal phase difference between the output of the delay itself and the 
local clock. Sufficient gain is provided to keep the phase tracking 
error less than 3 degrees. . 

After synchronization of the phases is acquired, change in the delay 
position is a continuous measure of phase change between the signal 
arriving from j and the 7 oscillator. A voltage proportional to this 
delay position eminates from the compensating buffer store (see box 
representing this store in Fig. 2). This voltage, designated 2;;, is trans- 
formed in a series of signal processing operations and used to control 
the i and j oscillator frequencies in such a way as to reduce the long 
term frequency differences and minimize the deviation of the compen- 
sating delays from their center positions. 

The control signal, x;;, may first be modified in a nonlinear circuit 
that emphasizes the large amplitudes in order to enhance the control 
influence of delays approaching overflow. Then it 1s fed separately to 
two summing circuits, one located at station 2 and the other at station 
j. These circuits generate weighted sums of all control voltages derived 
from links arriving at a station and subtract from the total a weighted 
sum of the control voltages from links departing the station. The con- 
trols returned from remote ends of the departing links are delayed 
in the simulator, as they would necessarily be in a real system. 

The output of the summing amplifier is filtered and applied to a 
saturating circuit that limits the maximum amplitude excursions. It 
is then used to adjust the oscillator frequency proportionately. 


III. DESIGN SPECIFICATION 


In designing simulators, it is unnecessary to copy specifications of 
realistic systems. Indeed, it is often advisable to use scaling factors 
that make the simulator convenient to operate. However, for this 
work, realistic specifications for a possible continental network are 
followed in order to gain practical experience in designing and operat- 
ing a network, and to discover any factors of practical importance 
that may have been overlooked in theoretical studies. 

An important component of a PCM center is the local clock which 
determines the timing of all switching actions at the center. In the 
simulator, no digital processing is performed, therefore, the clock is 
simply a 1 megahertz sinusoidal oscillator, one cycle of which rep- 
resents a frame of information in the digital system. A frequency con- 
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trol range of one part in 10° would be acceptable in practice. A suitable 
crystal tuned clock would typically have a long term stability of a 
little less than one part in 108%. This sets a limit of about 1 per cent on 
the realizable control accuracy. To maintain this accuracy, noise and 
zero drift introduced into the low-pass control signals must be less 
than 1 per cent of full amplitude. 

Transmission delays lie in the range 5 to 8 microseconds per mile, 
depending on the transmission medium. Thus, delays in a continental 
system may typically lie in the range 0 to 24 milliseconds. In typical 
conditions, such delays would vary by amounts up to 10 microseconds, 
largely because of temperature change. The synchronizing system 
compensating for this variation would have storage capacity in excess 
of the expected variation in order to have some control margin. 

Components of the simulator have properties similar to those dis- 
cussed above, but the range of parameter variation and disturbances 
have been increased in order to provide adequate margin for tests 
under extreme conditions. For measuring convenience, all sources of 
noise and drift in the system have been made much less than 1 per 
cent of signal, so that representative disturbances can be introduced 
in a controlled fashion. 

The servodrive for the compensating delay is designed so that its 
filtering action on the frequency control signal can be neglected. For 
this approximation to hold, the low-frequency linear response of the 
servo has been made equivalent to a delay of less than 2 milliseconds. 
Hysteresis and other tracking errors are less than 1 per cent of a cycle. 


IV. THE OSCILLATORS 


It is necessary to have very stable oscillators in order to keep 
spurious phase errors small. The 1 megahertz crystal oscillators pur- 
chased and modified for this application deviate, on average, less 
than one cycle of phase in ten minutes. This is the approximate dura- 
tion of an experiment. 

Hach oscillator is individually mounted in a temperature-stabilized 
enclosure, and is electrically screened from the other units. 

The oscillator frequencies are controlled linearly over a range +1 
cycle per second (that is, one part in 10°), by changing the voltage on 
a varactor diode that is loosely coupled to the crystal. The center 
frequency and voltage sensitivity are both adjusted by precise po- 
tentiometers connected as shown in Figure 3. Figure 4 gives an oscilla- 
tor control characteristic. . 

When the simulator is in use, the zero settings of the oscillators are 
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Fig. 3— Circuit for controlling the 1 MHz station oscillator. The oscillator is 
controlled linearly through a range of +1 hertz. 


checked before each experiment and their sensitivities are measured 
about once per week. Realignment, which seldom is necessary, is 
accomplished by observing Lissajou ellipses which portray the oscil- 
lator phases with respect to one another and to an atomic clock. 

Zero settings are adjusted to give stationary ellipses when the 
control voltage is zero. Next, the sensitivities are equalized by ad- 
justing the control voltage to give stationary patterns when all the 
oscillators are modulated simultaneously with a low-frequency sine 
wave. In this measurement one oscillator is used as a reference. Its 
absolute sensitivity is set using calibrated de inputs and timing beats 
with respect to the atomic clock. 
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Fig. 4— Normalized frequency response to a control signal for a typical station 
oscillator, 
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V. SIMULATION OF FIXED TRANSMISSION DELAY 


Fixed transmission path delays are simulated with commercially 
available magnetostrictive lines of 2, 6, and 8 milliseconds. Each 
link uses three of these in tandem. Thus, values of delay from 0 to 24 
milliseconds are available in increments of 2 milliseconds. Each delay 
line is accompanied by an amplifier that makes up the amplitude 
loss in the magnetostrictive line and its transducers. The net loss is 
about (40 + 37) dB where ¢ is the nominal delay in milliseconds. A 
representative delay unit is shown in Fig. 5. 

The gain of each of these amplifiers is tailored to match the as- 
sociated delay loss. A patching coaxial cable is used to couple the 
amplifier output to the next delay unit. A filter circuit is incorporated 
that reduces the effect of using cables of different lengths. 

To simulate the fixed delay in the return de contro] signal paths, 
lumped constant LC delays are used. Values of 2, 6, and 8 milliseconds 
are provided by z-sections selected in sets to equal the corresponding 
ac delay. 
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_ Fig. 5— Transmission delays are simulated with fixed magnetostrictive delay 
lines. Associated amplifier compensates for insertion loss. 
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VI. THE CONTINUOUSLY VARIABLE DELAYS 


Continuously variable magnetostrictive delay lines are used to simu- 
late the compensating buffer stores that would be used in a real PCM 
network. Similar delay lines are also used to introduce variation in 
the length of the simulated transmission links. 

Fig. 6 is a simplified diagram of the compensating buffer store. 
Each sonic delay wire is mounted horizontally and is paralleled by 
two resistive potentiometers. Potentiometer A is used to monitor the 
state of the delay. Potentiometer B is used to generate control signals. 
Terminal C is the output for the delayed ac signals. Correspondingly- 
labeled terminals are joined by a flexible cable. A sliding carriage 
holds the potentiometer wipers, a magnetostrictive receiving trans- 
ducer, and signal preamplifier. It is driven along a steel rail by a 400- 
cps motor through a gear train, windlass, and steel cable. 

The preamplifier, together with a main amplifier, compensates for 
approximately 66 dB of loss in the transducers. The loss in the sonic 
wire is insignificant; therefore, the output is virtually independent of 
delay position. The variable line has a larger transducer loss than a 
fixed line because of the inefficiency of the mobile receiver transducer. 

Torsional mode lines, rather than conventional longitudinal lines 
were selected for their greater delay per unit length, that is, approxi- 
mately 3.3 microsecond per centimeter. The maximum delay is a little 
more than 200 microseconds. Because of the need to resolve small 
fractions of a cycle (that is, 1 per cent), the electromechanical drive 
requires great precision. This is further complicated by the need, al- 
ready mentioned, for a short time constant in its linear response to 
slow disturbances. Design details of the servomechanism for the com- 
pensating delay shown in Fig. 6 are given in the Appendix. 

The “variable transmission delay” simulator has less stringent 
requirements. It is controlled by a conventional positioning servo- 
mechanism shown in Fig. 7. The motor drives these delays so that the 
potentiometer voltage at B approximately equals the applied voltage 
Vi. Six volts are needed to fully deflect the lines and the tracking error 
corresponds with an input of less than 50 mV. 

Fig. 8 shows a completed unit which is representative of the 24 
that are used for both transmission delay variations and for compen- 
sating buffer storage. 


VII. THE NONLINEAR FUNCTION 


Nonlinear shaping of the frequency control signal is conveniently 
accomplished by loading the potentiometer output at B with a 
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Fig. 6—The compensating buffer store is a servo-controlled, continuously 
variable magnetostrictive delay line. The maximum delay is slightly more than 
200 microseconds. 


variable resistor, as shown in Fig. 6. The control voltage is then 
given by 


ee (See 
dls) —e Ge V 
where 
po 24, 
b 


x is the fractional deviation of the delay from center, +V and —V 
are the voltages applied to the potentiometer, R, is the potentiometer 
resistance (5K), r is the load resistor which is variable from 220 
to 10K. 

Use of this potentiometer loading introduces contact noise into the 
control, but this has not been a serious problem during the 6 months 
the system has been in use. 
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Fig. 7— Variable transmission delay simulator. The variable delay line as- 


sembly is identical to the type used in the buffer store but differs in the control 
circuits. 


VIII. THE PHASE DETECTORS 


The phase detector circuit is intended to measure the small phase 
discrepancy between two supposedly locked sinusoids. For this pur- 
pose, a linear response is needed only over a small part of a cycle. 
However, it is important that the phase indication be independent 
of input amplitude. Otherwise, spurious amplitude changes, such as 
might be caused by delay line imperfections, would perturb the con- 
trol and could possibly cause instability. 

Fig. 9 shows the phase detector schematically. The two input sig- 





Fig. 8— One of 24 variable delay units. Twelve are used as buffer stores, the 
remainder serve as transmission delay simulators. 
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Fig. 9— Phase detector for measuring small phase discrepancies between two 
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nals are applied directly to individual waveform squaring circuits.® 
The phase difference between the resulting squarewaves is measured 
by determining the duration of their time overlap in an and gate. This 
gate is biased to give zero mean output when the two inputs are in 
quadrature, that is, have 50 per cent overlap, as seen in Fig. 10(a). This 
figure shows the net output which has been smoothed and fed through 
the de amplifier. The amplifier incorporates the filter needed to 
stabilize the servo, specifications of which are derived in the Ap- 
pendix. 

The fact that the system synchronizes with the signals in quadrature 
at each phase detector is inconsequential to the simulator experiments. 
A linear phase response is available for almost +90 degrees deviation 
from the quiescent phase. This is more than ample because 3 degrees 
phase displacement drives the motor to full speed. 

Fig. 10(b) shows that the phase indication is sensibly independent 
of input amplitudes within a large range. 
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Fig. 10—(a) Output of the phase detector as a function of phase difference 
between two input signals. (b) Phase errors resulting from amplitude variations 
in the local clock, v., and delayed signals, va, are negligible over the required 
operating range. 
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Fig. 11—Summing of weighted amplitudes. 


IX. THE OSCILLATOR CONTROL PATH 


The signals used to control the oscillator frequency are combined 
and shaped in a series of operations depicted here. The summing 
of weighted amplitudes is performed by conventional operational 
amplifier circuits, Fig. 11. Similar amplifiers are also used in the filter 
circuit, Fig. 12, and the amplitude saturator, Fig. 13. 

The properties of these circuits can be summarized as follows. The 
electrical gain in each limb of the summing amplifier is continuously 
variable over the range 1/11 to 11. The filter has a simple low-pass 
characteristic, 1/(1+pT), for signals below 100 cycles per second. 
The time constant JT is continuously variable from 0.2 to 200 seconds. 
The saturator clipping levels are continuously variable from 0.1 to 
10 volts and they automatically remain symmetrical about zero. 

The zero drift and low-frequency noise in these circuits is less than 
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Fig. 12— Low pass filters. 
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Fig. 13 — Amplitude saturator. 


50 mV at the oscillator control inputs and they are designed to ac- 
commodate 10 volt peak amplitudes. 


X. OPERATION OF THE SIMULATOR 


The simulator is arranged to have two active states: standby and 
run. In the standby condition, the external disturbances are discon- 
nected from all the variable delay simulators and from the oscillators. 
The compensating delay control servos are disconnected from the 
phase detector outputs and are connected instead as simple position 
servos. In this quiescent state the positions of all delays and the oscil- 
lator frequencies are determined by preset controls. 

To operate the system in the test mode two switches are pro- 
vided. One connects the externally-provided disturbance signals to 
the variable delays (Fig. 7) and to the oscillators. The other switch 
connects the compensating buffer stores (Fig. 6) in the normal run 
position. 


XI, REACTION TO FAILURES 


In a real PCM network it is imperative that a failure in any one 
section of the system not disable the entire network. One assurance 
against such a catastrophe is to provide an interlock that requires 
an output from both square-wave generators in the phase detector 
before the associated servodrive can take the delay from its quiescent 
setting. Then, if a transmission link fails or an oscillator signal is 
lost, all the associated compensating delays immediately go to their 
quiescent states and have no further influence on the frequency con- 
trol of the remaining operating system. 
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This facility is part of the simulator. It is used to simulate re- 
sponse in various fault conditions. 


AI. SUMMATION 


The simulator has been used to investigate a large number 
of networks and control parameters. It has performed reliably 
and functions close to theoretical prediction. Some specific results 
that have been obtained are reported in the article which follows 
this.t 

Figs. 14 and 15 give a general view of the complete simulator. 
Fig. 14 shows the array of 24 servodriven transports which serve as 
variable delay units and Fig. 15 shows the master control and record- 
ing equipment. 
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Fig. 14— Bays containing the 24 servo-controlled variable delay lines. 
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Fig. 15— Left to right: recording equipment, master control, and the fixed 
transmission delay simulators. 


APPENDIX 


Design of the Servodrive for the Buffer Stores. 


The gain and phase properties required in the servodrive loop are 
described here. Let Ax be the departure of the shuttle from perfect 
tracking. Then using capitals to denote Laplace transforms, X(s)/ 
AX (s) is the open loop gain, G, of the servo in Fig. 6 and the closed 
loop response will then be 


X(s) _ _ Gs) 
X(s) + AX(s) 14+ G(s) 
where 
_ _ GQ) 
OS) = Sa + sP) 
where 


Go = ab is the gain constant 
a is the volts per unit error in x out of the phase detector 
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bis the shuttle velocity per volt into the servo amplifier 
T, is the mechanical time constant of the motor and drive 
mechanism. 


The stabilizing filter has the form 


7 (1 + sT’) 
Ie)= (1 + sT.)(1 + sT)) 


Where the phase advancing time constant J’ is necessary to stabi- 
lize the control and the two retarding time constants ensure a falling 
gain characteristic at high frequencies. This falling gain is needed to 
reduce noise and to limit below 200 cps the band of signals fed to 
the 400 cps motor drive. 

The loop can be written 


1+G 
ae) 


1te(d49)4 5 & tT etD)+ © (T,T.+T2T3+T3T,) +s‘ {734 — 


and at low frequency as s tends to zero this resembles a simple delay 
G,* . It is concluded that a gain G) greater than 500 is needed to meet 
the delay specification of 2 milliseconds. But the condition for stability 
as given by Routh is 


TTLG,( 1" te 1) 
Gt @T TT, + TPT, +7, TD 


fed Ae -+ Pals + TT, 
Dy Tes 2s 








> ae ial 
Given 7, = 5 msec and selecting JT’ = 5 msec, T, = 1 msec, 73 = 
0.8 msec, the system is marginally stable with G) = 2250 and in practice 
was satisfactory for applications using gains up to 10°. 

Experimental measurements confirmed that the effective low fre- 
quency delay approximates 1/G,, it was also noted that tracking 
errors decreased approximately as 1/G) when the amplifier gain was 
varied. This tracking error is less than +0.003 centimeter on the 
magnetostrictive wire with a loop gain of 10°. 
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Organic Synchronization: Design of the 
Controls and Some Simulation Results 


By J. C. CANDY and M. KARNAUGH 
(Manuscript received October 24, 1967) 


Organic synchronization is a method for phase-locking the signals of an 
extensive digital communication network. Each clock in the network 1s 
made to depend on the phase drift of signals arriving and departing tts 
station. This work demonstrates the practicality of such schemes. A four 
station simulation operating in real time with realistic parameters for 
a transcontinental network is used to evaluate various types of linear 
and nonlinear controls and to study effects of changing clock frequency 
and transmission delays. Considerable attention is given to the analysis 
of linear organic systems in order to pave the way for reasonable choices 
of design parameters and to make the results more easily understandable. 
The experiments show that the systems are very stable and easy to implement. 
No difficulty was experienced in starting the systems or in modifying 
their structures and they were immune to large scale breakdown caused by 
local faults. 


I. INTRODUCTION 


A model for certain mutually-synchronized systems of clocks and 
transmission links has been described by M. Karnaugh.'. He calls 
systems that conform to this model “organic systems.” The work on 
organic synchronization has been motivated chiefly by a desire to 
synchronize the sampling and switching operations in a geographi- 
cally widespread pulse code modulation communication network. 

Broad sufficient conditions for the stability of nonlinear organic 
systems have not yet been mathematically established. Nevertheless, 
there is reason to believe that systems having readily achievable 
clock stabilities and transmission delays of terrestrial magnitude can 
be well behaved. Because of this, and because performance under a 
variety of starting conditions, parameter choices, and perturbations 
is of interest, an analog simulator for organic synchronization has 
been constructed and put to use. 
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Details of the simulator hardware are described in the preceding 
article.2 The present discussion will center on what has been learned 
from experimenting with the simulator. Preliminary attention is given 
to the analysis of linear organic systems; this aids the selection of 
reasonable control parameters for the design and study of some rep- 
resentative systems. 

The experiments have supported the conjecture that systems of 
continental, or even global, dimensions will be stable and easy to 
implement. No great difficulty in starting the systems or in modify- 
ing their structures has been encountered. 


Il. SYSTEM ORGANIZATION 


2.1 Discussion . 

The simulator contains four oscillators. Each of these corresponds 
to the local clock at a geographically distinct switching station. These 
stations may be interconnected through selected delay lines in any 
or all of the twelve possible directed paths. The sinusoids transmitted 
through these paths correspond in period to the data frames in a 
pulse code modulation communication system. 

Tach transmission path includes a large, fixed delay and a small, 
continuously variable, delay. The latter is used to simulate the slow 
variations in transmission delay which may occur in cables. In addi- 
tion, there is another continuously variable delay line at the receiv- 
ing terminal of each path. This simulates the buffer store which is 
needed to retime all data arriving at a switching station. The buffer 
store must synchronize the incoming data frames with the local 
switching actions. The latter are timed by the local clock. 

For example, consider the two stations illustrated in Fig. 1 which 
shows a signal transmitted from station j to station 7. The arriving 
signal is held in close phase agreement with the ith clock by the 
servoloop that automatically adjusts the buffer delay. Intuitively, 
we see that a constant frequency difference between the two clocks 
would drive the buffer at constant speed to one end, where synchron- 
ism would be lost. To prevent such failure, the clock frequencies are 
controlled by voltages derived from the position of the buffer, as 
shown. Gain factors a and b are placed in the frequency control paths. 

A complete network is composed of many links resembling the one 
we have described. When the frequency of the clock at each station 
depends only upon the states of the buffers in transmission links arriving 
at that station, 6 = 0, and the controls will be called “one-sided.” 
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Fig. 1— A transmission from station j to station 7. 


When the clocks are also controlled by the states of buffers in links 
leaving their stations, b ~ 0, the controls will be called ‘‘two-sided.”’ 
Two-sided controls require the transmission of narrowband control 
signals between nearest neighbors in the system. In the simulator, 
separate baseband delays are used for this purpose. 

The simulator incorporates filters, h, in the control paths to the 
clocks. They are used to explore the possibility of shaping the dynamic 
response of the network. Also included in the control paths is an ampli- 
tude limiter, p, which places a bound on the frequency deviations. 
Another nonlinearity, v, is placed in series with the buffer-position 
output, x. This causes the effective gain to vary with buffer position, 
so that the control influence of buffers near overflow may be made 
greater than those near their center position. 


2.2 Analysis 


Karnaugh has analyzed one-sided controls to determine the settling 
state after switching on from specified initial values. Now we shall 
examine changes in the settling conditions resulting from disturbances 
of transmission delays and clock frequencies, with two-sided controls. 
The model used by Karnaugh,’ with a very slight change in nota- 
tion, is illustrated in Fig. 2. Mathematically, 


f() = F; + | ncn > a, v(x; ;(t)] — py b; v(x, (¢ — 7.01} | (1) 


That is, the frequency of the ith clock equals its natural frequency 
F;, plus a control function of x, the buffer states at the 7** station and 
its nearest neighbors. 
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Fig. 2— Station 7 and its control with respect to station 7. 


The function p(-), which represents a limiter, is defined by 


p(t) = 2, if |x| SG 
= G, if «>G 
SOG ee ag (2) 


It is stipulated that the control filters, h, have unity gain for de: 


H,0) = [hl dt = 1, a Osea (3) 


where h;(t) is the impulse response and H;(0). is the de response. 
Mathematical convolution is indicated by ‘*.” 

The gain coefficients a;;, bi; are nonnegative values. They are both 
equal to zero when there is no transmission link to station 2 from sta- 
tion 7. When the link is present, a;; > 0 and, if and only if the controls 
are two-sided, bi; > 0. 

The function v(-) is monotonic and of odd symmetry. It is also 
assumed that 


v(+1) = +1. (4) 
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Simple examples of such functions are 


v,(2) = genth n= 0,1, Qrses (5) 


This is linear in the special case, n = 0. 

The instantaneous state of the buffer delay in the link to center i 
from station 7 is denoted by 2i;(t). Let the number of clock cycles 
stored in this buffer be y;;(t), and let its capacity be 2D,;. Then 


vist) = [yis() — D;\Di; . (6) 
This represents the fractional deviation of the buffer delay from its 


half capacity. For synchronism to exist. the buffers must lock the 
phase of the incoming signal to that of the local clock. 


z;(t) = D5 [pitt =) = OEP Ca (7) 
Where p;(¢) is the phase of the 7th clock, therefore 
pi(t) = f.(0) (8) 


The constant C,; is determined by the initial switch-on conditions.* 

7:4; 18 the transmission delay in the link to station 7 from station 7. 
7;; is the delay in sending control signals to station 7 from station j. 
It might or might not be approximately true that the two delays are 
equal when both are defined, but this has been true of the simulations 
that have been done. 

Now assume that the system has settled down to a common fre- 
quency f and that the phase differences remain finite. Thereafter 


pitt) = fi tr; (9) 


where 7; 1s independent of time. Also, because the buffers are quiescent 


z(t +7) = 2;,() forall 7 > 0 


and 
. RAL) ae Se 
so equation (1) becomes time independent, 
N 
f=F;,;+ p > {a:v(tii) — b:v(@;)}, (10) 
and (7) becomes 
as = Dol; — re — fra] + Cas, (11) 


which is also time independent. 
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For the purpose of investigating the changes in settling frequency 
of a system, F; and 7; will be considered variables, and such varia- 
tions are possible in the simulator. It is also possible to perturb the 
simulator by connecting or disconnecting transmission links, or by 
adding or removing an entire station. 


III. LINEAR SYSTEMS 


3.1 Analysis 


The linear subclass of systems is of very special importance. Study 
of these systems has added much to our knowledge of organic syn- 
chronization. 

One-sided linear systems have been very extensively studied. V. E. 
Benes, in unpublished work, was the first to establish a very inter- 
esting sufficient condition for stability of these systems, and to study 
their settling frequency. A. Gersho and B. J. Karafin®? have recently 
simplified the derivation and proof of these results. M. B. Brilliant‘ 
has shown that the requirement for network connectedness can be 
weakened, if master-slave relations between subsystems are per- 
mitted. He also has determined dynamic responses for some networks.® 
Karnaugh has formulated a model that makes explicit the influence of 
the starting conditions in the formulas for settling frequency.? 

Here, we shall extend the analysis to two-sided controls, consider- 
ing the network to be switched on and in a quiescent state and then 
determining the change of quiescent conditions resulting from certain 
parameter changes. The result is a linear expression for the small 
changes in frequency that result from small changes in delay and in 
natural frequencies of the clocks. 

The model defined in the previous section is linearized by removal 
of the limiter and by assuming v(x) linear for the small changes. 
Total derivatives of equations (10) and (11) show how the settling 
state defined by the frequency f and the buffer states x is affected 
by changes. Thus for small changes, A, from a given quiescent state 
we have 


N 
Af = AF; + pe {a;;0’(x;;) Av; — b;v'(2;,) Ax;,} (12) 
and 


Az;; = D;i[Ar; — Ar; — A(fr7;;))- (18) 
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Now use 
a3; = a,;D550' (2:3) 

and (14) 
By, = b;.D5iv' (eis) 


which are assumed constants for a given settling state. Also use 


a= Deu + By) (15) 
and 
r= Denes — Biers) (16) 
to get 7 
Af = AF; + >» (a;; + B;.) Ar; — g; Ar; — A(fr,). (17) 


Now measure phases with respect to the first station 


G=Di-D=ri—n, Ag = Ar, — An. (18) 
Then 


AF; — A(fri) (19) 


N 
Af + 2d {9 63; — (ai; + Bi} Ag; 
where 8,; is the Kronecker delta 


05; =] 63; = 0 for t x je 


I 


The solution of equation (19) by determinants, d, can be carried 
out to yield, after some simplification and index permutation, 


M= 


oy d; AF, B 4 (d; aig d; Bi) A(fri;) 
Af = a ae a eR ee ee (20) 


N N 
dX 4 dd; 
i=1 


t=1 


t 


This equation is not solved for Af explicitly because the term A (fzi;) 
depends on Af, however this result is a convenient one for practical 
use. Let us investigate its meaning. 

By definition, fz; is the phase delay in the transmission link con- 
necting center 7 to center 2. Therefore equation (20) expresses the 
change of settling frequency as a weighted sum of the changes in 
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clock natural frequencies, AF’, less a weighted sum of the change in 
phase delays during transmission, A (fz). 

The usefulness of the result derives from the fact that in most 
practical networks the phase delay variations are determined by 
changes of delay alone, variations resulting from changes of fre- 
quency, rAf, are usually negligible, that is 


A(fr:;) =f Ar; 


and this is true for the simulations. 
Returning to the general result, for small changes we can write 


A(fr:;) = f Ari; 1 T15 Af. (21) 
Then equation (20) becomes 


s d (AF, — f Ar;) 
1 ar (22) 


ye ad; ae T;) 


This is similar to the result given by Bene’. It demonstrates the 
surprising property that the sensitivity of the network to changes 
is small when the delays are large. However, the effect has only the- 
oretical interest because in practice 7; < 1. 

ach weighting coefficient, d; , is the cofactor of the element in the 
ith row, first column, of the matrix [T — XJ], defined by 


Nay = i; HO Bis 
Ds; = gs 9s;. 
It can be shown, following Brilliant,* that d; > 0 when 


N 
| S| >0 forall 7 #7 
n=1 qt 


and d; = 0 otherwise. The criterion for positive d; is, heuristically 
speaking, that a chain of transmission links of \ shall run from station 
7 to every other station. If the 7th station provides a master clock for the 
system, d; > 0 and d; = 0, j ¥ 7. The strongly connected systems 
of BeneS, on the other hand, have d; > 0 for all 7. The system cannot 
be synchronized unless at least one of these coefficients is positive. 

The solution, (22), of the linear equations depends upon their being 
nonsingular. This is equivalent to the nonvanishing of the denominator. 
This condition will be seen to hold in systems of practical interest. 
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3.2 The Expression for Setiling Frequency 


To give direction to the simulation experiments, we summarize the 
salient properties of equation (20) in a set of rules that can be readily 
confirmed experimentally without need for solving determinants. They 
are also important because they describe conditions that might be 
used in a real network. 


Rule 1 


In the case of two-sided controls that are “proportional” so that 
Kyja;; = K,@;; for allt # 7, where (K,K, --+ Ky) is a set of N positive 
numbers, then it can be shown that 


KG Kea Ks 
dy d, dy 
and that 
N 
>> K, AF; 
ON aes ear (23) 


N 
2d, Ki 


The settled frequency is a weighted average of the oscillator center 
frequencies and is independent of delay. 


Rule 2 


In the particular case of ‘‘balanced”’ two-sided control defined by 
a;; = B;; for all 7 # j, 


Af = D3 AF, . (24) 


=| 


The settled frequency of a balanced two-sided organic system is always 
the unweighted average of the clock center frequencies. 


Rule 3 


In the case of a reciprocal control defined by K;(a;; + 8;;) 
K;(a;; + 8;;) for allz ¥ j, it can be shown that 
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and that 


N N 


> K; AF; — D» {Kias; — KiBi;} AGfrs3) 


Aj = =) ______________. (25) 
K; 


Having reciprocal control requires that the controls also be linear, 
v(x) = x. One sided controls, 8 = 0, can be reciprocal if Kia = Kyjoj. 


Rule 4 

In the case of proportional control, as in rule 1, if a transmission 
link is symmetrical, 7; = 7; then, in response to a change Ary = 
Arji 


Az;; = Dj} Ar;jf and Az;; = Dj} Ar;.f 


and there is no change in the position of any other buffer in the 
network. 


Rule 6 


The net phase delay around a closed loop in a synchronized network 
is a constant integer. 
To prove this we notice that the net phase delay in link 7 is 


6:5 = Dil + 2:3) + rif. (26) 
Therefore 
Adi; = Di; Axi; + A(frii) 
using (13) 
Ad¢;; = Ar; = Ar; . (27) 
Then around a closed loop in the network the net phase delay will be 
Adi DOsn A Sy, Bone = Ary Arp Are = Arg 2s PDT; 
= 0. 


The definition of synchronization requires that the net phase delay 
be an integer. Notice that changes in the phase delays must be given 
a direction. That is, changes in delays introduced into signals flowing 
in one direction around the loop have the opposite sign from changes 
in delays introduced in signals flowing in the other direction. 
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IV. DYNAMIC RESPONSE AND STABILITY 


In the foregoing analysis, stability was tacitly assumed and the 
actual form of the dynamic response to the parameter changes was 
ignored. Here we shall explore conditions for stability of two-sided 
controls for a linear network having no limiter. The clock frequencies 
will be considered variables but for convenience we shall assume 
transmission delays are fixed. 

The transient response and stability of linear organic systems may 
be studied by means of the Laplace transforms of equations (1) and 

(7). Combining these equations and putting 


ri(t) = p:(t) — p-(0) (28) 
we obtain the result 


sk ,(s) = H,(s) > (a B;.)R,(s) 


— HRW Ys + B+ V9) QD) 


where capitals denote Laplace transforms, and 


Qi; = ae" (30) 
Bis = Bee (31) 
Bi = Been (32) 
V.) =*24 cw (33) 


and C,(s) represents initial conditions of the network. Following 
Gersho and Karafin,? we notice that equation (29) may be put in 
the form 


Bis) __ 
g:H ,(s) 


. Bus 


R.(s) = Bis) y fut be Gu TB Ry V.(9), (34) 


where 

gH (8) ; 
s+ HW) D (as + Bd 
Let M be the matrix whose (27)th element is 


Bs) = (35) 


M.,; = BG) artes (36) 
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and let Q be the diagonal matrix with 





Bs) = 
u= = 37) 
@ gH ;(s) ( 
The system is now seen to obey the vector equation 
R(s) = [I — M]"'QV(s). (38) 


When V(s) is specified, we can compute the solution R(s) from this 
relation. 

Arguments which differ from those of Gersho and Karafin only in 
minor detail can now be used to establish a sufficient condition for 
the stability of connected systems. The proof will not be repeated 
here. The resulting condition is that 

|Bs)| <1 for s=wV-14¥0 i=1,2,-::,N, (39) 
where w is real. 

B,(s) is independent of the system delays when the controls are 
one-sided. This is not the case for two-sided controls, as one may 
see from equations (32) and (35). However, the sufficient condition 
may be checked for any particular network and is easily satisfied 
in practice. 

M. B. Brilliant has simplified (39) for a simple but revealing case: 


H(s) = 1, gi: = 9, Tip HF = 7 for all t Fj. 


The sufficient condition for stability is then 


gr < 0.5. (40) 


Thus, the largest product of the delay and the gain can, if it is less 
than 0.5, guarantee stability of the unfiltered two-sided linear net- 
work. 


V. GAIN AND BUFFER DELAY 


It will be shown later that, for one-sided controls, limiting of the 
frequencies can cause synchronization failures when the transmis- 
sion delays are changing. Therefore, we are motivated to use small 
enough values of gains, g;, to avoid limiting. 

When there are no filters, or only single-pole filters, equations (1), 
(2), (4), and (6) show that limiting will not occur if 


> (ai; + b;) < G. (41) 


ORGANIC SYNCHRONIZATION 239 


Notice that G would be made equal to the maximum tolerable fre- 
quency deviation. For simplicity, suppose that 
g:=g9 for += 1, 2,---,N 
and 
D,;; = D forall (, 9). 
Then the inequality, (41), becomes 


gD <G. (42) 


Now consider the buffer size. The buffers deflect to compensate 
for delay changes and to correct oscillator drift. Let the greatest trans- 
mission delay variation from midrange be denoted by Ar and let the 
maximum error magnitude of the oscillator center frequency be denoted 
by AF so that 


AF = max | AF, |. 
Then we require 
D> FAr+ a (43) 


and we wish to keep the buffer sizes small by making g large. Then 
most of the buffer delay capacity is used to compensate for trans- 
mission delay variations 
F &r > . (44) 
Inequalities (44), (43), and (42) lead to 


F Ar D Par + 


Values of D that satisfy the last two inequalities in (45) will exist, 
provided that 


g< & = 
rs AF 
F Ar + — 
g 
whence, 
AF G — AF 


a OS =. (46) 
T T 
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Inequalities (46) are satisfied over a positive interval of values of 
g whenever 


G > 2 AF. (47) 


Finally, D must be chosen in the interval, 





Pane Spe = (48) 
g g 
For some specific examples, let us make the following choices: 
G=4AF 
* 2 AF 
F Ar 
D = #¥F Ar. 


Let us further assume the largest single transmission delay to be 
2 X 10~” seconds, which is of transcontinental magnitude, and let 
Ar = 2 X 10°° seconds. The latter is almost surely an overestimate 
for underground coaxial cable, but is reasonable if about 10 per cent 
of the cable is above ground. 

Two values of F will be used, corresponding roughly to voice and 
video sampling rates. For each of these, two values of AF will be used, 
corresponding roughly to the accuracies of simple crystal oscillators 
and atomic oscillators. Table I shows the resulting parameters. 

The greatest value of g encountered in Table I is 3X10. Let us 
examine the consequences of this with respect to the stability condi- 
tion, equation (40). The product 7g is less than 1210-* for all cases. 
Therefore in satisfying the stability condition 7g < 0.5 we have a 
factor of 400 to spare. 

From (35), we see that B;(s) is the transfer function for the ith 
clock with respect to equal perturbations of phase in all arriving 


TaBLE I—Parameters for Voice and Video Sampling 


F AF Ar G g D 
(Hz) (Hz) (seconds) (seconds!) (cycles) 
10! 3xX10-% 2X10 112X107 3X10 3.5107 
108 107 2X 10-8 4X 10-5 1074 3.5 X 1071 

6x10° 118X107 2X10° 72109 3X10% 21X10 
6 X 108 6X10-- 2X10 24 X 107 10-4 2.1 X 10 
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signals. In the absence of delays, 


_ Gg: H(s) 
2S s + g:H,(s) . 
Therefore, if H;(s) is either constant or is monotone decreasing with 
little phase shift in 


s=avV-—l, Ofw89:, 


the response radian bandwidth is no more than g;. When H; = 1, 
the response has the simple time constant, g;’. 

This enables us to interpret 7g as the ratio of the delay between 
stations to their time constants of response. When this ratio is very 
small, the transmission delays have negligible effect on the dynamic 
response of the network. This is another strong reason for using rela- 
tively low gains. 

Another useful approximation that applies to most practical net- 
works concerns the phase delay during transmission fr;;. The change 
of phase delay is 


A(fr:;) = f Aris + 75; Af. (21) 


We have seen that delays change about 0.1 per cent* while clock 
frequencies change about 10-’, thus the second term in the expression 
is negligible in most cases and this is true for the simulations. Equa- 
tion (20) then becomes 


> d (AF; — f Ar;) 
Af 2 =—__-______ (49) 


N 

de ds 

t=1 
VI. A TWO STATION NETWORK 


6.1 Response of the System to Small Change 


Many revealing properties of the organic synchronizing scheme 
can be illustrated with two stations connected in a loop as in Fig. 3. 
Such a system is characterized by the following equations 


Af & Gi APs + gc AP; (asc: — B:sBii)f Ares + tis) 


50 
gi + 9; gi + 9: ( ) 


* For short delays the fractional change can be much larger. 
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Fig. 3— Two coupled stations, 7 and j. 








Aitsy SS ALF, — Fe + (ais + Bit Altar + 70 
gi T Qi (51) 
~N Dj = 
Ax; = Pees [AIP; — Fi] + (@is + Bid f Aas + 7ia)). 


Equation (50) was derived from (49), and equations (51) follow from 
(13), (18), (19), and (49). 

Our simulation of this network has a nominal 1 megahertz center 
frequency with +1 hertz control range. The gains (g; + g;) are in 
the range 0.1 to 0.001 sec”*. Transmission delay around the loop can 
be preset in the range 0 to 0.1 second and be varied continuously by 
+100 microseconds in each link. 

In all setups the settling states agree well with prediction. They con- 
firm that fixed delay has negligible effect (« one percent) on buffer 
and frequency deflections. Therefore the approximation in deriving 
(49) from (20) is justified. 

For illustration, Table II gives some typical results. The data have 
been normalized to represent unit amplitude disturbances. The in- 
cremental delay changes, Av are expressed as a fraction z of the asso- 
ciated buffer capacity such that 

Ar;,F 
Sean Ye (52) 

In the simulations the gain delay product 7g is less than 10-*. This 
should guarantee stability with reasonable valued filters, and indi- 
cates that fixed delays are too small to have appreciable effect on 
transients. Observations of responses after switch-on and of subse- 
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Taste II—Settling States of Two Stations 


Gain Setting Disturbance Result 
Control Type ais Biz ay = Byi Fs AF; Azji Aziz Axjt Af 
Balanced 1 1 1 I 1 0 —1/4 41/4 1/2 
1 1 1 1 0 —1 1/2 1/2 0 
2 2 2 2 1 0 —1/8 +1/8 1/2 
2 2 2 2 0 = 1/2 1/2 O 
2 2 1 1 1 0 —1/6 —1/6 1/2 
2 2 1 41 0 = 1/3 2/3 O 
Proportioned 2 1 1 2 1 0 —1/6 1/6 1/3 
2 1 1 2 0 —1 1/2 1/2 0 
1 2 2 1 1 0 —1/6 1/6 2/3 
1 2 2 1 0 —1 1/2 1/2 0 
One Sided 1 0 1 0 0 —1 1/2 1/2 1/2 
2 0 2 0 1 0 —1/4 1/4 1/2 
2 0 2 0 0 —l 1/2 1/2 1 
2 0 1 0 0 —1 1/3 2/3 2/3 


Unit gain = 10-? sec! 

Unit frequency = 1 Hz 

Unit delay = 10? usec 

Nominal clock frequency = 1 MHz 


quent disturbances confirm this conclusion. Indeed, delay had insig- 
nificant effect when various one- and two-pole low pass filters were 
included in the control loops. The filters investigated had cut-off 
frequencies from 0.01 to 50 hertz and q-factors up to 10. 

When fixed delays are neglected, the response of the linear system 
to small change can be expressed as , 


AG F(s) — F;(s) — {a;.(s) + B;.(s) [7.;(8) + 79 ]F 
De B+ 7.0) + 9.0) we 
f(s) & a;;(s)X.;(s) — B;;(s)X;,(s) (54) 
where 
a; ;(8) = a;;H,(s), B,;(s) = 8.;H ;(s) 
and 
gi(s) = g:H;(s). 


When no filters are used these responses have a simple time constant 
1/(g: + 9;). Fig. 4 gives some typical response curves for a system 
with no filters and using various gain values. Including delays up to 
0.1 second had no noticeable effect on these curves. 
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fe) eo 50 75 100 125 150 
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Fig. 4 — Response of buffer a; to a step decrease in delay Ari;. The controls 
are balanced, and the net gain is such, that (g:+g;)~ is 12 seconds for curve a, 
24 for b, 48 for c, 108 for d, and 216 for e. 


Fig. 5 shows responses with similar one-pole low-pass filters in 
both oscillator circuits. Notice that use of filter time constants near 
1/g speeds the response. We do not anticipate that filters will play an 
important role in the:operation of the networks. We have seen that 
they are not needed to stabilize the controls, nor will they be needed 
to speed responses because most changes will be thermally induced, 
and will therefore occur very slowly. Some filtering may be needed 





























ie) 50 100 150 200 250 300 350 
SECONDS 


Fig. 5— Response of buffer x:; to a step decrease in 71; with various filters 
1/(1+s7). The controls are balanced and the net gain (g:+g:) = 1/24 sec.*. The 
filter time constants 7 are 0 seconds for curve a, 10 for b, 20 for c and 40 for d. 
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to smooth the output of phase detectors and to reduce noise. For this, 
time constants not greater than 1/g should be adequate. 


6.2 When Control Signals Limit 


A saturating limiter p is placed in the control path to the oscillators 
in order to place a bound on the frequency deviations. To demon- 
strate its effect, Fig. 6 shows a buffer response to a ramp change of 
delay using various gain values in a symmetrical one-sided control 
system with no filters 


aii(s) = a;:(s) = 9 Bii(s) = B,:(s) = 0. 

At the start of the experiment the network is quiescent with the 
two clock natural-frequencies offset from one another and the two 
buffers defected by amounts that provide sufficient control voltage 
to align the running frequencies. 


Az;; = AP i D;; . 
gi 

This initial buffer-displacement decreases proportionately with 
increased gain. For this purpose a large gain is attractive. However, 
when the transmission delay varies, limiting occurs if the gain is 
large, as shown in curve c of Fig. 6, where there is loss of control 
and ultimately loss of synchronization. Clearly, saturation must be 
avoided, that is, equation (41) satisfied, when one-sided controls are 
used. 






CORRECTION FOR LIMITING 


INITIAL FREQUENCY<< 
OFFSET 
-0.25 





7 
- 





BUFFER 
a OVERFLOWS 

















ie) 25 50 75 100 125 
SECONDS 


Fig. 6 — Response of zi; to similar ramp delay increases Arij and Arj:; in a 
system having initial frequency discrepancy. The controls are one-sided and the 
ee such that g:' and g;“ are both equal to 50 seconds for curve a, 25 for 

and 5 for c. 
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With two-sided “proportioned” control, changes of delay do not 
cause change of settling frequency. See rule 1. Any frequency changes 
that might occur during the dynamic response will usually be small 
because large parameter changes will occur slowly in a practical net- 
work. Thus with “proportioned” control, delay changes will not cause 
limiting, the signal in the limiters being determined only by the dif- 
ferences of oscillator center frequencies. 

Larger gains can be used with two-sided contro] systems than with 
one-sided. However, in practice equation (41) should always be 
obeyed so that systems continue to operate reliably if faults inter- 
rupt the distant control paths. The reasons for avoiding limiting are 
equally valid for more complex networks, and this is not a hard 
restriction on design. It has already been seen in Table I how limit- 
ing can be avoided while using gains that satisfy the main re- 
quirements of a practical network. 


6.3 Nonlinear Control 


Use of nonlinearity in the control loop has been proposed’ as a 
means for exaggerating the correcting influence of buffer stores that 
are near overflow, at the expense of those nearer center. For this purpose 
the nonlinear circuit v(x) is included in Fig. 2. It makes the incremental 
gains a function of buffer position, as illustrated by equation (14). 
We shall examine separately the nonlinear response to incremental 
delay and incremental frequency change as given by equations (50) 
and (51). 


-0.25 











LINEAR 












NONLINEAR 








100 125 


50 75 
SECONDS 


_ Fig. 7— The response of both buffers to a ramp increase in delay Ari; with 
linear and nonlinear controls. Controls are balanced with net g = 0.04 second* 
and p= 1.2. 
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When 2z;; is nearer zero than 2;; the nonlinearity can be selected 
to make | v’(;;) | very small compared with | v’(z;;) |. Then (a;; + B:;) « 
(a;; + 8;;). Therefore x;; will approximately track the delay change 
and x;,; will have little change. 

To illustrate this, Fig. 7 shows the response of both buffers to a 
steadily increasing delay commencing with z;; = 0 and 2;; & 0.5. 
The nonlinearity successfully makes the buffer with the most reserve 
compensate for most delay change. In this and subsequent experiments 


v(x) = Ger ; p> 1. (55) 


A useful method for demonstrating the response of nonlinear control 
is a locus on a graph of 2;; plotted against x;; . Fig. 8 shows such loci 
for various starting conditions, and Fig. 9 is an enlargement of a section 
of Fig. 8 repeated for different values of nonlinearity. These results 
show that the nonlinearity can prevent overflow in some circumstances. 

Next, consider the response to frequency change with constant delays. 
Equation (51) shows that when the deflections are expressed in cyles 


Lij —_ 


_-\atl=1 OVERFLOW 


BOUNDARY 


NONLINEAR 


FINISH 





fe) 


Fig. 8— Loci of xj; against x:; for a slow ramp change of delay Arij in a two- 
station network; 1: balanced linear controls gi = gj = 0.04 second’; n: bal- 
anced nonlinear controls w = 1.2 gi = g; = 0.04 second”. 


248 THE BELL SYSTEM TECHNICAL JOURNAL, FEBRUARY 1968 


ctij—_ 
OVERFLOW BOUNDARY 


FINISH 
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Fig. 9— An enlargement of the upper left quadrant of Fig. 8 for various values 
of nonlinearity. 


they are the same at both stations. We need to minimize their value 
when a buffer is near overflow. Consider a worst case with z;; = 0 
and x;; & 1. Then there is improvement over the linear control if 
the net gain (g; ++ g;) exceeds the linear gain. That is, if 


[v’(0) — 1](a;; + Bi) + [v’') — 1](ai; + Bis) > 0. 


For example, when (a;; + B:;) & (ai; + 6;;) there is improvement 
if 1 < » < 8 using nonlinearities described by equation (48) or if 
n > 2 using nonlinearities described by equation (7). An illustration 
is given in Fig. 10. 

The degree of nonlinearity used will be limited by the need to maintain 
some minimum gain for controls near center. 


VII. A CHAIN OF FOUR STATIONS 


The simulation is extended by adding two more stations to form 
the chain shown in Fig. 11. In designing controls for this network it 
seems sensible to make controls associated with each link reciprocal 
and symmetric, that is, 


(ai; + Bi:) = (ai; + 8s;) forall ¢ ¥ 3. 
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Fig. 10— Loci of 2;; against x; for a step change in frequency AF; in a two- 
Gatien network with balanced nonlinear control uw = 1.2. 


Then with the buffers at center the settling frequency is equally 
sensitive to frequency drift in each oscillator. With such control the 
net gain, g, of the relay stations j and k exceeds that of the end sta- 
tions 2 and l. The relay stations therefore have greater risk of limiting, 
but this risk is lessened by the larger possibility of averaging and 
canceling of disturbances. 

Some typical settling states for this chain are given in Table III. 
The top of Fig. 12 shows how a sudden change of station 2 frequency 
causes transients through the system. The bottom shows the corre- 
sponding response with nonlinearities added. Notice the slowing down 
of the response and the increased deflections particularly at remote 
stations. These effects are a consequence of the low control gains 
when buffers are near center. 


AT; L 
eC ME * ) 


Fig. 11 — The four station chain. 
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Fig. 12 — Buffer responses to a step change of frequency AF, in a four station 
chain. Top: linear balanced controls a = 8 = 0.02 second. Bottom: nonlinear 
balanced controls a = 8 = 0.02 second” » = 1.2. 


However, we have said that frequency drift is small in practice, 
more significant are effects of delay change. These are demonstrated 
in Fig. 18 which shows loci of x; against zy, for a delay change Arj. 
For simplicity, buffer x;; starts with the same content as x. The con- 
tent of x; and 2, are apparent from the graph; all other buffers start 
at center. Curve a is the response with two-sided balanced control. 
For slow change, it is a vertical line whose shape is independent of 
position on the graph or of nonlinearity value. Curve b is the response 
with one-sided linear control. The others are for nonlinear one-sided 
controls. 


VIII. VARIOUS CONNECTIONS OF FOUR STATIONS 


Figs. 14, 15, and 16 show four stations connected as a loop, star, 
and complete network, respectively. The settling states of these three 
networks are illustrated in Tables IV, V, and VI. These results were 
virtually independent of fixed delay values and confirm the predic- 
tion of equation (49). 
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Fig. 13 — Loci of x;; against 2, for a slow ramp change of delay Arv;; in a four 
station chain. 


Fig. 14— Ring network. 
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Fig. 15 — Star network. 
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Fig. 16 — Complete network. 


Fig. 17 gives some typical transient responses of the loop connec- 
tion for both linear and nonlinear controls. Fig. 18 demonstrates the 
effect of nonlinearity on loci of x; to x; for various starting positions. 
A corresponding graph for the fully connected network is shown in 
Fig. 19. 

In studying the responses of all these networks no noticeable effects 
were observed when delays up to 0.1 second were included in signal 
and control paths. The responses were always stable even with single 
pole filters in the control loop. 

Obtaining proper statistical data on the response of the synchron- 
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_ Fig. 17 — Buffer responses to a step change of frequency. AF; in a four station 
ring. Top: balanced linear controls g = 0.012 second. Bottom: nonlinear con- 
trols » = 1.2 g = 0.012 second”. 
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TaBLE IV — Settling States of a Four Station Ring 
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ized network to disturbances is an extensive project because of the 
large number of interacting parameters available and because of the 
slow response of the system to change. To short cut this work we 
have examined tendencies in the response in some extreme conditions. 

The fully connected network was disturbed by driving the delays 
and the oscillators individually with triangular waves having unre- 
lated frequencies in the range 0.01 to 0.001 Hz. The maximum deflec- 
tion of the twelve buffers and their likelihood of overflow under 
various control settings were observed. 

The following conclusions were made from these experiments. 


(2) The chance of buffers overflowing is decreased by use of in- 
creased control gain, provided limiting is avoided. Increasing the 
gain past the limiting point increased the chance of overflow. 

(uw) The chance of overflow increases when limiting levels are 
decreased. 

(2) The chance of either limiting or overflow occurring is less 
with balanced control than with one-sided control. 

(tv) Use of low pass filters with long time constants increases the 
chance of overflow. | 

(v) Use of nonlinearities (u & 1.2) in a balanced control system 
reduces the chance of overflow. 

(vz) Use of nonlinearity with single ended control increases the 
chance of overflow. 

(viz) Use of nonlinearity increases the chance of overflow when 
disturbances change at rates comparable with the control bandwidth. 


IX. RESPONSES IN FAILURE 


If a single transmission link is severed or a signal is lost for any rea- 
son the servo in Fig. 1 normally drives its buffer to an extremity and 
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TABLE V — Settling States of a Four Station Star 
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thus deflects the frequency of the whole network. In the simulator this 
is prevented by auxiliary equipment that takes the buffer to center 
in case of signal loss. Meanwhile, the remainder of the network can 
function correctly. 

Another interesting fault occurs when a buffer is driven to over- 
flow by the normal controls. In the simulator it rests at full deflec- 
tion until the phase drift brings another frame within half a cycle 
of synchronism with the local clock, the system locks to that frame. 

Little difficulty has been experienced in switching stations in and 
out of an active network. The switching transient can be reduced by 
adjusting the station frequency so that its buffers are near center. 


X. CONCLUSION AND RECOMMENDATION 


There seems to be no good reason for using nonreciprocal con- 
trols, moreover using linear-symmetric controls (aj + Bu) = (az + 
By) makes the network equally sensitive to frequency drift in each 
oscillator (rule 3). Then if the oscillators have similar properties this 
minimizes the frequency displacement caused by oscillator drift. 








0 


Fig. 18 — Responses of buffers to a ramp change in delay Ar;; in a four station 
ring: (a) linear balanced control; (b) linear one-sided control. All other curves 
for nonlinear balanced control » = 12 starting with vj, = 22), v143 = vj: and all 
other buffers at center. 
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Fig. 19— Responses of buffers to a ramp change of delay Ar;, for the com- 
plete network with the same controls as in Fig. 18. 


There is additional advantage in also balancing controls (aij = Bi). 
Then slow delay change will not affect the system frequency and the 
chance of buffers overflowing is reduced. The disadvantage of using 
two-sided, and hence also of balanced controls, is the need for control 
paths between centers. These paths require little bandwidth (< 1 Hz) 
but must pass de; they could be included in the framing and signal- 
ing codes of a PCM system. 

A commercially useful network will almost certainly contain both 
sending and receiving links in each path used, and these links will 
be similar in opposite directions. For example, they will expand and 
contract together so that A7;; = Ar. Rule 4 shows that if the controls 
in such a network are proportional then each delay change is cor- 
rected directly by the buffer in its link and no disturbance propagates 
through the network. 

If the buffers in such a symmetric path start at corresponding posi- 
tions they will approximately track one another, deviating only to 
correct small frequency drift and the misbalance of the links. There 
will be little use for nonlinear shaping to equalize their positions, 
besides, nonlinearities lower sensitivity when buffers are near center 
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and increases the practical difficulty of keeping controls proportional. 

Use of symmetric links and linear balanced controls provide for 
near optimum design of buffers. Each one corrects only the delay 
change of its own link and the small frequency drifts. They can be 
designed independently of the network. 
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AC Focusing of an 
Electron Microscope Objective Lens 


By R. D. HEIDENREICH 
(Manuscript received November 29, 1967) 


Control of phase contrast in electron microscope images is presently 
achieved by fine vernier control of the de objective lens current. An 
incremental change 8J in lens current produces a corresponding in- 
crement Sf in focal length and a focal separation in object space of 
6L,. L, is the distance from lens to conjugate object plane so that AL, 
is the axial separation of conjugate object plane and scattering speci- 
men. The defocus phase at the image plane of a ray scattered at angle 
B relative to the axial ray is 


K 
x= Blane (1) 
with |K| = 27/,. The phase contrast G in the image goes as 


Gxcie an Lat ALB? (2) 


ignoring spherical aberration, astigmatism, etc., for this purpose. Pre- 
cise control of phase contrast for a spacing “a’ = AB requires that 
the focal separation increments 8L, be as small as 100 A and prefer- 
ably less. 2 

For unsaturated magnetic lenses the change in focal length with 
lens current is taken as® 


if = —20,,% (3) 


with C',, the chromatic coefficient. The focal increment 6Ly) ~ 6éf in 
this case. For saturated lenses 6Z, and 6f may not be equal so that it is 
important to calibrate 6L, in terms of the objective vernier control 
using a Bragg beam from a crystal. 

As a result of re-examining the effect of drift and ripple in the high 
tension and in the objective lens current, the possibility of achieving 
a selective or “tuned” focus for certain desired object spacings suggested 
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itself. If the lens current has a harmonic ripple sin wt, then the focal 
separation AL, is likewise a periodic function of time and so is the phase 
contrast (2). The response of the regulated de supply to an ac voltage 
applied across the lens through a condenser apparently results in a 
time-dependent focal separation that behaves like 


AL (t) = —e(1 + sin wi) (4) 


with g, the amplitude in angstroms proportional to the applied ac 
voltage. The time average phase contrast obtained by using (4) in 
(2) is found to be 


G),~ J {E ob) sin “ (AL, — ¢%.)8°. (5) 


If yg, is zero, the Bessel function is unity and there is just the usual 
phase contrast. J, has appreciable values at its subsidiary maxima or 
side bands being —0.4, +0.3, —0.28, etc., occurring at (K/2) 98? = 
3.8, 6.9, 10.1, etc., respectively. Since it is very difficult to reduce the 
ripple in a power supply below about 1 mv it is seen that for certain 
values of residual ripple it is feasible that the Bessel function in (5) 
could be zero for certain object distances. The injection of an ampli- 
tude ¢, of the right amount could give “side band contrast” in J, for 
a particular distance “a” in the object. 

Fig. 1 is an eek of the effectiveness of this type of focusing. 
These photographs are from a focal series taken of a graphitized 
carbon black particle and exhibit the 3.4 A (002) spacing of graphite. 
Growth defects in the form of dislocations and curvature of the planes 
are evident. The ac injection increments in Fig. 1 are 1.5 mv, 2.5 mv 
and 5 mv, with the corresponding amplitudes 9, roughly 300, 500, 
and 1000 A, respectively. Both the Bessel function and the de effect 
of y. are operative in the contrast. The micrographs were taken with 
a Siemens Eleniskop I modified for short focal length (2.4 mv at 100 
KY). 

Since the ripple in the high tension is an uncertainty the absolute 
values of g, to achieve maximum side band contrast are uncertain. 
A continuously variable ac potential is thus desirable with the lens 
voltage being monitored by an oscilloscope. 

The possible advantage of “side band” or “microphase” focusing 
lies in obtaining useful contrast at spacings a < 2 A. To fall well within 
the principle maximum of J,, the residual ripple for a 1. 5 A spacing 
must be below 1 part in 10°. The injection of 9, & 75 A, 125 A, or 
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Fig. 1— Electron micrographs showing the 3.4 A (002) spacing in graphitized 
carbon. black. By injecting 1.5, 2.5, and 5 mv of 60 cycle ac into the lens coil, the 
focal series shown by a, b, and c, respectively, was obtained. The center picture 
is very near optimum and shows growth defects in the lattice, whereas a is un- 
derfocused and ¢ slightly overfocused for the 3.4 A spacing. 
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225 A ac should “tune” the respective side bands for the 1.5 A spacing. 
The (004) period of 1.7 A in graphite has been imaged in this way 
using carefully aligned, azzal illuminations. 
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Air-Insulated Beam-Lead Crossovers 
for Integrated Circuits 


By M. P. LEPSELTER 
(Manuscript received December 4, 1967) 


Air-insulated crossovers for integrated circuit interconnections have 
been developed which feature low capacitance and dual dielectrics 
for high yield and reliability. This type of crossover is fabricated 
using beam lead technology and is applicable to integrated circuit 
chips as well as substrates. 

The cross section drawing in Fig. 1 shows a crossover with dual 
dielectric insulation consisting of a 5 micron air gap in series with 
a 2000 A layer of ZrOo. Fig. 2 is a plan view photograph of a typical 
pattern. The dual dielectric crossover will not fail under mechanical 
load even if the top conductor is deflected until it touches the bottom 
level, since the ZrO. has a breakdown strength by itself of approxi- 
mately 100 volts. (See Fig. 3.) Any pinholes in the solid dielectric 
will not lead to short circuits either, as they would if the top level 
metallizing were applied by deposition. When the external force is 
released, the top beam lead will elastically return to its original posi- 
tion, and the air gap will once again sustain over 200 volts. 

The fabrication procedure is as follows. A Ti-Pt lower level contact 
pattern is formed according to methods previously described.? A layer 
of Zr-Cu is then deposited, and feed-through holes etched where con- 
tact to the lower level is desired. The top layer, a gold beam lead, is 
electroformed using the copper as a base, and the copper is etched 
away. Thermal oxidation of the Zr layer to ZrOz completes the struc- 
ture: two levels of metallization separated by a ZrOs-air dielectric. 

The use of a thick air gap as one of two dielectrics between two 
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Fig. 1— Cross section of crossover. 
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Fig. 2— Plan view of crossover pattern on integrated circuit chip (narrowest 
line is 0.4 mil). 


crossing conductors leads to a system where pinhole shorts are nearly 
impossible, thermal stresses are minimal, and parasitic capacitance 
is low (~0.001 pf/mil?). Conventional crossovers consisting of evap- 
orated metals over and under a deposited dielectric (for example, 
$i02) are prone to pinhole shorts, conductor opens, and high thermal 
stresses owing to the mismatches involved. Diffused crossunders are 
limited to use in chips only, and use large areas of silicon because of 
their high sheet resistances. 

The beam lead crossover system has been successfully applied to 
integrated circuit chips as well as complex substrates (see Fig. 4). 
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Fig. 3— Cross section of loaded crossover. 
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Fig. 4— Photograph of integrated circuit memory substrate. 


In addition to the high initial yields, aging of hundreds of crossovers 
in room air at 350°C for over 1000 hours has shown no apparent 
change, nor is any envisaged. The metallurgical system has been 
stringently tested on beam lead transistors and found to be very 
rugged. It is therefore believed that these crossovers will be generally 
useful in many integrated circuit applications. 
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